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ABSTRACT 

We present peculiar velocities for 84 clusters of galaxies in two large volumes at 
distances between 6000 and 15000 kms~^ in the directions of Hercules-Corona Bore- 
alis and Perseus-Pisces-Cetus. These velocities are based on Fundamental Plane (FP) 
distance estimates for early-type galaxies in each cluster. We fit the FP using a maxi- 
mum likelihood algorithm which accounts for both selection effects and measurement 
errors, and yields FP parameters with smaller bias and variance than other fitting 
procedures. We obtain a best-fit FP with coefficients consistent with the best exist- 
ing determinations. We measure the bulk motions of the sample volumes using the 
50 clusters with the best-determined peculiar velocities. We find the bulk motions in 
both regions are small, and consistent with zero at about the 5% level. The EFAR 
results are in agreement with the small bulk motions found by Dale et al. (1999) on 
similar scales, but are inconsistent with pure dipolc motions having the large ampli- 
tudes found by Lauer & Postman (1994) and Hudson et al. (1999). The alignment 
of the EFAR sample with the Lauer & Postman dipole produces a strong rejection 
of a large-amplitude bulk motion in that direction, but the rejection of the Hudson 
et al. result is less certain because their dipole lies at a large angle to the main axis 
of the EFAR sample. We employ a window function covariance analysis to make a 
detailed comparison of the EFAR peculiar velocities with the predictions of standard 
cosmological models. We find the bulk motion of our sample is consistent with most 
cosmological models that approximately reproduce the shape and normalisation of the 
observed galaxy power spectrum. We conclude that existing measurements of large- 
scale bulk motions provide no significant evidence against standard models for the 
formation of structure. 

Key words: galaxies: clustering — galaxies: distances and redshifts — galaxies: 
elliptical and lenticular, cD — galaxies: fundamental parameters — cosmology: large 
scale structure of universe 



1 INTRODUCTION 

This paper reports the main results of the EFAR project, 
which has measured the peculiar motions of clusters of 
galaxies in two large volumes at distances between 6000 and 
15000 kms"\ The project was initiated in the wake of early 
studies of peculiar motions which found large-scale coherent 
flows over significant volumes of the local universe (Dressier 
et al. 1987, Lynden-Bell et al. 1988). The primary goal of 



the EFAR project was to test whether such large coherent 
motions were to be found outside the local volume within 
6000 km s"^ In the following years, the velocity field within 
6000 km s"^ has been mapped by several methods and in 
increasing detail so that today there is fair agreement on 
the main features of the motions (recent results are given in 
Giovanelli et al. 1998a, b, Dekel et al. 1999, Courteau et al. 
2000, Riess 2000, da Costa et al. 2000, Wegner et al. 2000 
and Tonry et al. 2000; see also the review by Dekel 2000). 
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The bulk velocity within this volume, and its convergence 
towards the frame of reference defined by the Cosmic Mi- 
crowave Background (CMB), appear to be consistent with 
the broad range of currently-acceptable cosmological models 
(Dekel 2000, Hudson et al. 2000). 

However on larger scales there have been measurements 
of bulk motions which, at face value, appear much greater 
than any acceptable model would predict. The first of these 
was the measurement by Lauer & Postman (1994), using 
brightest cluster galaxies, of a bulk motion of ~700kms~^ 
towards (Z, 6)«(340°,+50°) for a complete sample of AbcU 
clusters out to 15000 kms^^. More recently, large motions 
have been also be obtained for two smaller samples of clus- 
ters at similar distances, for which peculiar velocities have 
been measured by the more precise Fundamental Plane and 
TuUy-Fishor estimators: Hudson ct al. (1999) find a mo- 
tion of 630±200kms"^ towards («, 6) = (260°,-l°) for the 
SMAC sample of 56 clusters at a mean distance of ^8000; 
Willick (1999) finds a motion of 720±280kms"^ towards 
(Z,6)=(272°,-hl0°) for the LPIOK sample of 15 clusters at 
very similar distances. Those two motions are in good agree- 
ment with each other, but are nearly orthogonal to the 
Lauer & Postman motion (though similar in amplitude). 
In contrast, the other extant study of peculiar motions on 
scales greater than 6000 km s"'^, the SCH TuUy-Fisher sur- 
vey of Dale ct al. (1999a), finds a bulk flow of less than 
200kms~^ for a sample 52 Abell clusters with a mean dis- 
tance of ~11000kms~^ 

At these scales the robust prediction of most cosmo- 
logical models is that the bulk motion should be less than 
300 km s"'^ with about 95% confidence. It is therefore of 
great interest to determine whether there really are largo co- 
herent motions on scales of ^10000 kms^'^. The EFAR pecu- 
liar motion survey probes the velocity field in the Hercules- 
Corona Borealis and Perseus-Pisces-Cetus regions, which 
are almost diametrically opposed on the sky and lie close 
to the axis of the bulk motion found by Lauer & Post- 
man. With 84 clusters in these two regions extending out 
to ~15000kms~^, the EFAR sample is well-suited to test- 
ing for this particular bulk motion. Conversely, however, it 
is not well-suited to testing for a bulk motion in the di- 
rection found for the SMAC and LPIOK samples, which is 
almost orthogonal to the major axis of the EFAR sample. 
The main goal of this paper is to determine the peculiar 
motions of the EFAR clusters and the consistency of the 
bulk motion of the sample with both theory and other bulk 
motion measurements on similar scales. 

The structure of this paper is as follows: In §2 wo sum- 
marise the main features of the data presented in Papers I- 
IV of this series. In §3 we describe the maximum likelihood 
gaussian algorithm developed in Paper IV, which is used 
to determine the parameters of the Fundamental Plane and 
obtain the distances and peculiar velocities for the clusters. 
In §4 we derive the best-fitting Fundamental Plane and crit- 
ically examine the random and systematic uncertainties in 
the fitted parameters. In §5 we derive the distances and pe- 
culiar velocities for the clusters, testing them for possible 
systematic biases and comparing them to the peculiar ve- 
locities obtained by other authors for the same clusters. In 
§6 we determine the bulk motion of the sample and compare 
it, using a variety of methods, to the results of other studies 



and to theoretical expectations. Our conclusions arc given 
in §7. 

We use l/o=50kms~^ Mpc and go=0.5 unless otherwise 
specified. All redshifts and peculiar velocities are given in the 
CMB frame of reference. 



2 THE EFAR SAMPLE AND DATA 

Earlier papers in this series have described in detail the 
selection of the clusters and galaxies in the EFAR sample 
(Wegner et al. 1996, Paper I), the spectroscopic data (Weg- 
ner ct al. 1999, Paper II; Colless ct al. 1999, Paper V), the 
photoelectric and CCD photometry (Saglia et al. 1997a, Pa- 
per HI; Colless et al. 1993) and the photometric fitting pro- 
cedures (Saglia et al. 1997b, Paper IV; Saglia et al. 1993). 
In this section we briefly summarise the main properties of 
the EFAR database. 

The clusters of galaxies in the EFAR sample are selected 
in two large, distant (i.e. non-local) volumes: Hercules- 
Corona Borealis (HCB, 40 clusters, including Coma) and 
Perseus-Pisces-Cetus (PPC, 45 clusters). These regions were 
chosen because they contain two of the richest supcrcluster 
complexes (excluding the Great Attractor/Shapley super- 
cluster region) within 20000 km s~^. The clusters come from 
the ACO catalogue (Abell et al. 1989), the fist of Jack- 
son (1982) and from scans of Sky Survey prints by the 
authors. The nominal redsliift range spanned by the clus- 
ters is 6000kms"Vc«<15000kms"^ The distribution of 
the EFAR clusters on the sky is shown in Figure 2 of Par- 
pcr I; their distribution with respect to the major superclus- 
ter complexes is shown in Figure 3 of Paper I. 

Galaxies were selected in each cluster for their appar- 
ently elliptical morphology on Sky Survey prints, and for 
large apparent diameter. The total sample includes 736 
early-type galaxies in the 85 clusters. Apparent diameters 
were measured visually for all early-type galaxies in the clus- 
ter fields. The range in apparent visual diameter (Dw) is 
from about 10 arcsec to over 60 arcsec. The sample selec- 
tion function is defined in terms of these visual diameters; 
in total, Dw was measured for 2185 early-type galaxies in 
the cluster fields. Selection functions are determined sepa- 
rately for each cluster, and are approximated by error func- 
tions in log Dw ■ The mean value of the visual diameter is 
(log _Di4/)=1.3 (i.e. 20 arcsec), and the dispersion in log Dw 
is 0.3 dex (see Paper I). 

We obtained 1319 spectra for 714 of the galaxies in our 
sample, measuring redshifts, velocity dispersions and the 
Mg6 and Mg2 Lick linestrength indices (Paper II). There 
are one or more repeat observations for 45% of the sample. 
The measurements from different observing runs are cali- 
brated to a common zeropoint or scale before being com- 
bined, yielding a total of 706 redshifts, 676 velocity dis- 
persions, 676 Mg6 linostrongths and 582 Mg2 linestrengths. 
The median estimated errors in the combined measure- 
ments arc Ac2=20kms"\ Ao-/cr=9.1%, AMg6/Mg6=7.2% 
and AMg2=0.015 mag. Comparison of our measurements 
with published datasets shows no systematic errors in the 
redshifts or velocity dispersions and only small zeropoint 
corrections to bring our hnestrengths onto the standard Lick 
system. 

We have assigned sample galaxies to our target clus- 
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ters (or to fore/background clusters) by examining both the 
line-of-sight velocity distributions and the projected distri- 
butions on the sky (Paper II). The velocity distributions 
were based on EFAR and ZCAT (Huchra et al. 1992) red- 
shifts for galaxies within 3 h^Q Mpc of the cluster centres. 
These samples were also used to derive mean redshifts and 
velocity dispersions for the clusters. The original selection 
was effective in choosing cluster members, with 88% of the 
galaxies with redshifts being members of sample clusters and 
only 12% lying in fore/background clusters or the field. The 
median number of galaxies per cluster is 6. 

We obtained R-band CCD photometry for 776 galax- 
ies (Paper III) and B and R photoelectric photometry for 
352 galaxies (CoUess et al. 1993). Comparison of the CCD 
and photoelectric photometry shows that we have achieved 
a common zero-point to better than 1%, and a photometric 
precision of better than 0.03 mag per measurement. Circu- 
larised galaxy light profiles were fitted with seeing-convolved 
models having both an i?^''* bulge and an exponential disk 
(Paper IV). We find that only 14% of the galaxies in our 
sample are well fitted by pure R^^* bulges and only about 1% 
by pure exponential disks, with most of the sample requir- 
ing both components to achieve a good fit. From these fits 
we derive total R-band magnitudes rriT, Dn diameters (at 
20.5 mag arcsec"^), half-luminosity radii Re, and average ef- 
fective surface brightnesses (SBe), for 762 galaxies. The total 
R magnitudes span the range mT=10. 6-16.0 ((mT) = 13.85), 
the diameters span Dn— 4.8-90 arcsec {{D„}=2Q arcsec), 
and the effective radii Re span 1.6-71 arcsec {{Re)=6-9 arc- 
sec). For 90% of our sample the precision of the total mag- 
nitudes and half-luminosity radii is better than 0.15 mag 
and 25% respectively. The errors on the combined quantity 
FP = log Re—0.3{SBe) which enters the Fundamental Plane 
equation are always smaller than 0.03 dex. The visual selec- 
tion diameters Dw correlate well with the D„ diameters (or, 
equivalently, with the Fundamental Plane quantity FP). 

The morphological type classifications of the galaxies, 
based on all the information available to us, reveal that 31% 
of the sample objects, visually selected from photographic 
images to be of early type, are in fact spiral or barred galax- 
ies. The 69% of galaxies classified as early-type can be sub- 
divided into 8% cD galaxies, 12% E galaxies (best-fit by a 
pure R^^^ profile), and 48% E/SO galaxies (best-fit by a disk 
plus bulge model). 

All the EFAR project d ata is available from NA SA's 
Astrophysical Data Centre (littp://adc. gsfc.nasa.gov) and 
the Centre de Donnees astronomiques de Strasbourg 
(http//:cdsweb.u-strasbg.fr). A summary table with all the 
main parameters for every galaxy in the EFAR sample is 
available at these locations as J /MNKAS / vol / page. The 
contents of the summary table are described here in Table |l| 



3 MAXIMUM LIKELIHOOD GAUSSIAN 
METHOD 

We use a maximum likelihood gaussian algorithm for fit- 
ting the FP and determining relative distances and peculiar 
velocities. This algorithm, which is described in detail in 
Paper VI, was developed in order to deal with the general 
deficiencies of previous approaches and with some specific 
problems posed by the selection effects and measurement 



errors in the EFAR sample. Previous methods for fitting 
the FP using forms of multi-linear regression have not fully 
dealt with the intrinsic distribution of galaxies in size, ve- 
locity dispersion and surface brightness, nor with the simul- 
taneous presence of measurement errors with a wide range 
of values in all of these quantities. The maximum likelihood 
gaussian algorithm properly accounts for all these factors, 
and also handles complex selection effects in a straightfor- 
ward way. The selection criteria for the EFAR sample are 
well-determined, and involve both the original sample selec- 
tion based on galaxy size and a posteriori limits imposed 
on both galaxy size and velocity dispersion. A specific prob- 
lem with the data is that the velocity dispersion measure- 
ments include a significant fraction of cases where the er- 
rors, though themselves well-determined, are large relative 
to the actual value. There is also the fact that the numbers 
of galaxies observed per cluster are relatively small, so a 
method is required which is both efficient and robust against 
outliers (either unusual galaxies or errors in the data). The 
extensive simulations carried out in Paper VI demonstrate 
that the maximum likelihood gaussian method is superior to 
any of the classical linear regression approaches, minimising 
both the bias and the variance of the fitted parameters, and 
performing well in recovering the FP parameters and pecu- 
liar velocities when presented with simulations of the EFAR 
dataset. 

The maximum likelihood gaussian method assumes that 
each galaxy i is drawn from an underlying gaussian distri- 
bution in the three-dimensional FP-space (r = log Re, s = 
\oga,u = (SBe)). We also assume that this underlying dis- 
tribution is the same for each cluster j, apart from a shift 
Sj in the distance-dependent quantity r resulting from the 
cluster's peculiar motion. We want to determine the mean 
values (r , s, u) and the variance matrix V which characterise 
the galaxy distribution, along with the shifts Sj due to the 
clusters' peculiar velocities. We do this by maximising the 
likelihood of the observed galaxy data over these parame- 
ters, while properly accounting for all the various selection 
effects. 

The probability density for the ith galaxy, in terms of 



S,Ui — U), IS 



^^^'^ " (2^)3/211/ + 1 1/2/, " 



(1) 



where V is the variance matrix of the underlying distribu- 
tion and Ei is the error matrix of the measured quantities. 
The errors are convolved with the intrinsic dispersion of the 
galaxy distribution to give the observed distribution of the 
data. The exclusion function 6{y) = Y[ where 8{y) = 1 
if J/ > and otherwise, accounts for parts of FP-space that 
are inaccessible because of selection effects. For simplicity, 
we assume that these selection effects apply to linear com- 
binations of the variables, described by the matrix A. The 
normalisation factor fi is such that J P{x) dPx = 1, and ac- 
counts for the selection effects described the exclusion func- 
tion O. The likelihood of the observed sample is 



(2) 



where S{xi) is the selection function giving the probability of 
selecting a galaxy with parameters Xi. In order to correct for 
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1. Description of EFAR summary data table. 



Column 


Code 


Description [units] 


1 


GIN 


Galaxy Identification Number 


2 


CID 


Cluster Identification (see Paper I) 


3 


CAN 


Cluster Assignment Number (see Paper II) 


4 


Clus 


Cluster Name (corresponds to CID) 


5 


Gal 


Galaxy Name 


6 


RAh 


Right Ascension (J2000) [hours] 


7 


RAm 


Right Ascension (J2000) [minutes] 


8 


RAs 


Right Ascension (J2000) [seconds] 


9 


Deed 


Declination (J2000) [degrees] 


10 


Deem 


Declination (J2000) [minutes] 


11 


Decs 


Declination (J2000) [seconds] 


12 


I 


Galactic longitude [degrees] 


13 


b 


Galactic latitude [degrees] 


14 


Tvne 


Morphological type 


15 


Dn 


Diameter enclosing a mean R-band SB of 20.5 mag arcsec"^ [arcsec] 


16 




Error in Dn [arcsec] 


17 


D„(20) 


Diameter enclosing a mean R-band SB of 20.0 mag arcsec"^ [arcsec] 


18 


D„{19.25) 


Diameter enclosing a mean R-band SB of 19.25 mag arcsec""^ [arcsec] 


19 




Half-luminosity radius in the R-band [arcsec] 


20 


-R^fkDc) 


Half- luminosity radius in the R-band [kpc, i^o— 50, ^o~0.5] 


21 


SBe 


R-band surface brightness at Re [j^i^g arcsec"^] 


22 


(5SBe 


Photometric zero-point error on SBe [niag arcsec"^] 


23 


(SBe) 


Mean R-band surface brightness inside Re [niag arcsec"^] 


24 


<5(SBe> 


Photometric zero-point error on (SBe) [mag arcsec"^] 


25 




Total apparent R magnitude [mag] 


26 


(5m 7^ 


Photometric zero-point error on mjp [mag] 


27 




Bulge half-luminosity radius in the R-band [arcsec] 


28 




Bulge R-band surface brightness at RqB [mag arcsec"^] 


29 


h 


Disk scale-length in the R-band [arcsec] 


30 




Disk central surface brightness in the R-band [mag arcsec ~'^] 


31 


' ''/ "-eis 


Ratio of bulge half-luminosity radius to disk scale-length 


32 


D/B 


Disk-to- bulge ratio (ratio of luminosity in disk to luminosity in bulge) 


33 


Fit 


Type of fit (B=bulge, D=disk, BD=bulge-|-disk; other, see Paper III) 


34 


p 


Quality of the photometric zero-point (P=0 good, P=l badj see Paper III) 


35 


Q 


Global quality of the photometric fit (l=best, 2=fair, 3=poor; see Paper HI) 


36 


B-R 


B — R colour [mag] 


37 


5{B-R) 


Error in B — R colour [mag] 


38 




Mean ellipticity inside Re 


39 


An 


Reddening in the R-band 


40 




Cluster mean redshift [kms~^] 


41 




Error in cZei [kms~^] 


42 


cz 


Galaxy redshift [kms~^] 


43 


Scz 


Error in cz [kms~^] 


44 


a 


Central velocity dispersion of galaxy [kms~^] 


45 


Sa 


Error in cr [kms~^] 


46 


Mgb 


Mgb Lick linestrength index [A] 


47 


SMgb 


Error in Mg6 [A] 


48 


Mg2 


Mg2 Lick linestrength index [mag] 


49 


<5Mg2 


Error in Mg2 [mag] 


50 


Qs 


Spectral quality (A=best, E=worst; see Paper II) 


51 


a/e 


Absorption/emission flag 


52 


log Dw 


Logarithm of the D\y diameter [arcsec] 


53 


S{Dw) 


Selection probability computed using D\y (see §H) 


54 


log Dw{D„) 


Logarithm of Dw computed from Dn (see §^ [arcsecl 


55 




Selection probability computed from DiY{Dn) (see gSl) 



The s ummary table is availahl e as 3 /MNHAS/ vol /page from NASA's Astrophysical Data Centre 
(ADC. httD: //adc.gsfc. nasa. gov ) and from the Centre de Donnees astronomiqucs de Strasbourg (CDS, 



http:/ /cdsweb.u-strasbg.fr). 
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the selection function, each object in the sample is included 
in the likelihood product as if it were 1/S{xi) objects. 

The error matrix can be computed from the estimated 
errors {Svi, 5si, SFPi, SZPi), where 5FP is the error in the 
combined quantity FP — r otu (with a ~ 0.3) and SZP is 
the photometric zeropoint error. In terms of these quantities, 
the error matrix for galaxy i is 



( 



\ 





f)Srj-SFPl 

a(l + Q,^) 













5u^ 



(3) 



Note that Ssi combines the estimated random errors in the 
velocity dispersion measurements and the correlated errors 
between galaxies introduced by the uncertainties in calibrat- 
ing dispersions obtained in different observing runs to a com- 
mon system (see Paper II). Likewise, Sui is given by the 
quadrature sum of the error on the effective surface bright- 
ness (from the fit to the galaxy's surface brightness distri- 
bution) and the photometric zeropoint error (see Paper III) . 

l)5fP? + + a^)5r^ 



(a^ 



+ SZPt 



(4) 



For the EFAR sample, the selection function depends 
on galaxy diameter and varies from cluster to cluster (see 
Paper I). For galaxy i, a member of cluster j, the selection 
probability is 



S, = S{logDv 



- erf 



iog7:ii. 



log^'Si-, 



.(5) 



The selection function for cluster j is characterised by -DxVj , 
the size at which the selection probability is 0.5, and by Swj, 
the width of the cutoff in the selection function. For early- 
type galaxies, the visually estimated diameter Dwi corre- 
lates with the measured diameter D„i according to the rela- 
tion log D„i — 0.80 log Dwi+0.26, with a scatter of 0.09 dex 
in log D„i (see Paper III). Because the visual diameters given 
in Paper I are individually uncertain, in computing selection 
probabilities we actually use an estimate of Dwi obtained by 
inverting this relation and inserting the accurately measured 
value of D„i. 

In order to avoid biasing the FP fits and the estimated 
peculiar velocities, it would be desirable to sample the same 
part of the FP galaxy distribution in all clusters. However, 
because the clusters are at different redshifts, the approx- 
imately constant apparent diameter selection limit corre- 
sponds to actual diameter selection limits D'^^ for the clus- 
ters that vary by about a factor of 2-3 (the approximate 
range of cluster redshifts; see Paper 1). We can limit this 
redshift-dependent sampling bias by excluding the smaller 
galaxies, which are only sampled in the nearer clusters. 
Guided by the simulations of Paper VI, we choose a selec- 
tion limit 73vi/cut = 12.6 kpc. This choice balances the reduced 
bicis of a higher Dwcut against the larger sample size of a 
lower Dwcut (95% of galaxies in the EFAR sample have 
£^vi/i> 12.6 kpc). Because of the good correlation between 
Dwi and the combined quantity FP — r — 0.3u (see Pa- 
per III), this cut in Dwi corresponds to an approximate 
selection limit FP^ut ~ 0.78 log Dwcut - 6.14 ^ -5.28. 

Another selection limit is due to the difficulty of measur- 
ing velocity dispersions smaller than a spectrograph's instru- 
mental resolution. For the spectrograph setups we used, only 



velocity dispersions greater than about lOOkms"^ could be 
reliably measured (see Paper II). We therefore impose a 
limit Scut=2, excluding galaxies with cr<100kms~^. The 
overall exclusion function for the EFAR sample is thus 
e = e(s-Scut)e(FP-FPcut). 

The mean of the distribution, (r, s, u), the variance 
matrix V , and the shifts Sj, are all determined by minimising 
— ln£, which for the EFAR sample is given by 



-\nC = 2_]sr'[0.5f,^(V + S0^« + 0.51n|V + £;,l+ln/.](6) 

FP>FJ=cut 

(where the constant term 1.51n(27r) has been dropped). The 
normalisation fi is obtained by integrating the gaussian dis- 
tribution over the accessible volume defined by s>Scut and 
FP>FPcut- The minimisation is performed using the sim- 
plex algorithm (Press et al., 1986). 

The FP is defined as the plane r = as+bu+c that passes 
through (r, s, u) and whose normal is the eigenvector of V 
with the smallest eigenvalue. For convenience, we define the 
second axis of the galaxy distribution to be the unit vector 
within the FP that has zero coefficient for s (in fact, this 
turns out to be a reasonable approximation to one of the 
remaining eigenvectors of V^). The three unit vectors giving 
the axes of the galaxy distribution can then be written in 
terms of the FP constants as 



vi = f — as — bu 
V2 = f + u/b 



(7) 



V3 = -f/b-{l + b)s/{ab)+u 



where f, s, and u are the unit vectors in the directions of 
the FP-space axes. The eigenvalues of V give the dispersions 
(Ti, (T2 and (T3 of the galaxy distribution in the directions of 
the eigenvectors; the smallest eigenvalue, ai, is the intrinsic 
dispersion of the galaxies about the FP. 

The final step of the process is to recover each clus- 
ter's distance and peculiar velocity. The mean galaxy size, 
r = log Re , provides a standard scale which we can use to de- 
termine relative distances and peculiar velocities. The offset 
5j between the true mean galaxy size, log Re, and the mean 
galaxy size observed for cluster j, log 7?e— is a measure of 
the ratio of the true angular diameter distance of a cluster, 
Dj, to the angular diameter distance corresponding to its 
redshift, D(zj): 



dex(log Re 



D{z,) dex(logPe 



10*^^ . 



(8) 



The relation between angular diameter distance and redshift 
(Weinberg 1972) is given by 

cz 1+ z + VI + 2i7o2 



D(z) 



(9) 



Ha{l + zY 1 + qoz + VI + 2(7oz ' 

We assume _ffo=50kms'^ Mpc, 50=0.5, and compute all 
redshifts and peculiar velocities in the CMB frame of ref- 
erence. The peculiar velocity of the cluster, Vj, is then ob- 
tained as 

_ czj-cz{Dj) 

1 + .(/),) ' ^'"^ 

where z{Dj) is the redshift corresponding to the true dis- 
tance Dj through the inverse of equation (^. Note that 
we are not using the low-redshift approximation V — cz ~ 
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HqD — cz{l — 10*), which leads to small but systematic er- 
rors in the peculiar velocities (e.g., at C2=15000 kms~^, the 
approximation leads to a systematic peculiar velocity error 
of about -4%). 

These distances and peculiar velocities are relative, be- 
cause the standard scale is determined by assuming that the 
distance (or, equivalently, peculiar velocity) of some stan- 
dard cluster (or set of clusters) is known. Distances and 
peculiar velocities are therefore in fact relative to the true 
distance and peculiar velocity of this standard. 



4 THE FUNDAMENTAL PLANE 

4.1 Best-fit solution and random errors 

We determine the parameters of the Fundamental Plane and 
the cluster peculiar velocities in a two-step process. We first 
fit the Fundamental Plane using only those clusters with 
6 or more suitable galaxies having reliable dispersions, ef- 
fective radii and mean surface brightnesses (the criteria are 
given below). We exclude clusters with fewer members be- 
cause the simulations of Paper VI show that including less 
well-sampled clusters increases the variance on the FP pa- 
rameters. We then determine peculiar velocities for all the 
clusters in a second step, where we fix the FP parameters 
at the values determined in the first step. This procedure 
results in more accurate and precise peculiar velocities than 
a simultaneous global solution for the FP parameters and 
the peculiar velocities. 

In order to be included in the fit a galaxy had to satisfy 
the following criteria: (1) good quality photometric fit (Q=l 
or Q=2; see Paper III); (2) a > 100 km s"^ and 51ogcr < 0.5 
(see Paper II); (3) a selection diameter Dw > 12.6 kpc and a 
selection probability >0.1. The first criterion excludes galax- 
ies with unreliable structural and photometric parameters 
(see Paper III); the second excludes galaxies with disper- 
sions less that the typical instrumental resolution or which 
have very large uncertainties; the third ensures that the clus- 
ters have uniform selection criteria and that no individual 
galaxy enters with a very high weight. No galaxy is excluded 
on the basis of its morphological type. There were 31 clus- 
ters in the sample with 6 or more galaxies satisfying these 
criteria. 

As well as these a priori criteria, we also rejected a fur- 
ther 8 galaxies on the basis that they lie outside both the 
3-CT ellipse of the galaxy distribution in the _FP-log a plane 
when the FP fit is obtained using all the galaxies in these 
31 clusters meeting the selection criteria (including them- 
selves) , and outside the 5-cr ellipse of the galaxy distribution 
when the FP fit is obtained excluding them. These galax- 
ies are listed in Table ^, which gives their galaxy ID num- 
ber (GIN), their cluster assignment number (CAN), their 
EFAR name, their morphological type and, where appro- 
priate, their NGC/IC numbers. The reasons why these 8 
galaxies are poorly fitted by the FP distribution that satis- 
factorily represents the other 255 galaxies fulfilling the se- 
lection criteria are not apparent. Although three are spirals, 
the other five include two ellipticals, two E/SOs and a cD. 
Three are members of A2151, including the cD NGC 6041. 
Two of these galaxies (GINs 45 and 370) are in clusters with 
data for 6 members; these two clusters (A160 and A1983) 



Table 2. Galaxies excluded from the Fundamental Plane fits. 



GIN 


GAN 


Name 


Type 


NGG/IC 


(i) Galaxies 


in clusters w 


ith >6 


members 


45 


7 


A160 C 


E/SO 




167 


21 


A400 H 


E/SO 




370 


43 


A1983 2 


S 




396 


46 


J16-W B 


s 




456 


53 


A2147 D 


E 




495 


58 


A2151 A 


cD 


N6041 


500 


58 


A2151 F 


S 


11185 


501 


58 


A2151 G 


E 


11193 


(ii) Galaxies 


in clusters with <6 


members 


355 


42 


J14-1 D 


S 




489 


57 


J18 C 


E 




552 


63 


A2162-S G 


E/SO 





therefore drop out of the sample of clusters to which we 
fit the FP. Also listed in Table ^ are another 3 galaxies in 
clusters with fewer than 6 members that are excluded from 
further analysis because they lie outside the 5-cr ellipse of 
the best-fitting galaxy distribution. 

The final sample of 29 clusters used to fit the FP pa- 
rameters is listed in Table ^, which gives the cluster as- 
signment number (CAN), the cluster name, the mean he- 
liocentric redshift and the number of galaxies that enter 
the FP fit. Of these 29 clusters, 12 are in HCB and 17 in 
PPC. They span the redshift range 6942 km s"^ (Coma) to 
20400 kms"^ (A419), though most are in the range 9000- 
15000 kms^^. However, they have similar selection diame- 
ters D'^, with minimum values of the Dw diameter in the 
range log Dvk (kpc) = 1.0-1.3. The Coma cluster sample is 
supplemented with the data of Miiller (1997; see also Miiller 
et al. 1998, 1999), which were obtained using essentially the 
same methodology. Midler's photometric data have been ad- 
justed by adding 0.04 mag in order to bring them into agree- 
ment with the EFAR data for galaxies in common. 

In fitting the FP we assume _ffo=50kms~^ Mpc~^ and 
go=0.5. We fix the zeropoint of the FP by forcing the mean 
of the FP shifts of the 29 clusters to be zero — i.e. we fix 
log Re by requiring ^5j=0. This results in a peculiar ve- 
locity for Coma of only — 29kms~^, so our FP zeropoint is 
essentially identical to that obtained by setting the peculiar 
velocity of Coma to be zero, as is often done. The effec- 
tive radii and mean surface brightnesses used were the total 
Re and (SBe) (rather than the bulge-only R^b and (SBes)) 
given in Paper III. In applying absorption corrections (taken 
to be 2.6i?B-v/4.0) we have adopted the mean of the ab- 
sorption corrections derived from Burstein & Heiles (1982, 
1984; BH) and Schlegel et al. (1998, SFD; with Eb-v off- 
set by —0.02 mag, the mean offset from BH given by SFD). 
The above assumptions and cluster/galaxy selection criteria 
yield our best fit to the FP. This best fit is given as case 1 
in Table ^ which lists the number of clusters and galaxies 
in the fit, the FP coefficients a, b and c, and the means and 
dispersions describing the galaxy distribution: logi?e, log a, 
(SBe), ai, (J2 and us. The table also explores the effects of 
the various assumptions and selection criteria, giving the FP 
fits obtained for a wide range of alternative cases. 

Case 1 is our best-fit solution. The EFAR FP, based 
on 29 clusters and 255 galaxies, has a=1.223±0.087, 
fe=0.336±0.013 and c=-8.66±0.33. The intrinsic scatter 
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Table 4. The parameters of the Fundamental Plane derived 



Case 






a 


b 


c 


logi?e 


logo- 


1 


29 


255 


1.223 


0.3358 


—8.664 


0.7704 


2.304 


2 


29 


271 


1.286 


0.3439 


—8.975 


0.7621 


2.298 


3 


29 


261 


1.201 


0.3265 


—8.430 


0.7840 


2.306 


4 


29 


255 


1.232 


0.3373 


—8.721 


0.7686 


2.300 


5 


29 


255 


1.183 


0.3292 


—8.422 


0.7961 


2.315 


6 


29 


255 


1.220 


0.3349 


—8.639 


0.7739 


2.306 


7 


29 


235 


1.235 


0.3357 


—8.690 


0.7750 


2.300 


8 


29 


255 


1.082 


0.3221 


—8.062 


0.8159 


2.329 


9 


29 


275 


1.132 


0.3224 


—8.184 


0.7773 


2.297 


10 


29 


244 


1.300 


0.3388 


—8.906 


0.7446 


2.292 


11 


29 


222 


1.247 


0.3303 


—8.607 


0.7369 


2.265 


12 


29 


255 


1.077 


0.3014 


—7.665 


0.7511 


2.310 


13 


29 


256 


1.207 


0.3359 


—8.625 


0.7725 


2.299 


14 


29 


258 


1.204 


0.3472 


—8.846 


0.6745 


2.206 


15 


29 


241 


1.080 


0.3239 


—8.081 


0.8099 


2.315 


16 


29 


255 


1.221 


0.3309 


—8.553 


0.8302 


2.331 


17 


29 


255 


1.223 


0.3345 


—8.629 


0.7895 


2.305 


18 


29 


255 


1.215 


0.3342 


—8.628 


0.7700 


2.307 


19 


29 


255 


1.227 


0.3359 


—8.648 


0.7990 


2.302 


20 


29 


255 


1.226 


0.3359 


—8.704 


0.7418 


2.303 


21 


29 


255 


1.227 


0.3374 


—8.707 


0.7735 


2.304 


22 


29 


255 


1.247 


0.3341 


—8.694 


0.7721 


2.301 


23 


66 


397 


1.206 


0.3274 


—8.452 


0.8021 


2.307 


24 


52 


355 


1.208 


0.3272 


—8.460 


0.7969 


2.306 


25 


39 


304 


1.244 


0.3265 


—8.531 


0.7927 


2.306 


26 


31 


265 


1.228 


0.3329 


—8.616 


0.7839 


2.305 


27 


16 


173 


1.109 


0.3432 


-8.525 


0.7487 


2.299 


28 


7 


99 


0.992 


0.3526 


-8.425 


0.7652 


2.326 


29 


29 


222 


1.330 


0.3351 


-8.904 


0.7776 


2.320 


30 


66 


348 


1.284 


0.3327 


-8.737 


0.8186 


2.330 


31 


52 


310 


1.293 


0.3323 


-8.756 


0.8097 


2.330 


32 


39 


267 


1.352 


0.3291 


-8.835 


0.7966 


2.323 


33 


31 


232 


1.333 


0.3300 


-8.804 


0.8020 


2.322 


34 


29 


223 


1.147 


0.3198 


-8.174 


0.7558 


2.300 


35 


29 


199 


1.241 


0.3125 


-8.250 


0.7568 


2.319 



various cases. 



(SBe> 




Ol 


0-2 


03 


Notes 


19.71 


0, 


.0638 


1.995 


0.6103 


st&rid&rd fit 


19.72 


0, 


.0688 


1.958 


0.6201 


includes Q— 3 photometry 


19.74 


0. 


.0671 


2.057 


6202 




19.73 


0. 


.0642 


1.992 


0.6138 




19.69 


0. 


.0632 


2.019 


0.5901 




19.71 


0. 


.0638 


1.996 


0.6043 


iicjpGi fi4.SFr)-l-n "^fiRH rnrrprtinnq 


19.74 


0. 


.0642 


2.014 


0.6161 


fivi^liiriPQ nl/~ici'jT^n 1 


19.74 


0. 


.0612 


2.057 


0.5394 


Tin JlTir r'lit 1 c; n nr^ 1 1 pH 


19.73 


0. 


.0675 


2.122 


0.6827 


TiTir J. — f\ 'ili"r\r' 

^ W cut — ^•*-> i^y^ 


19.69 


0. 


.0637 


2.040 


n 6220 


^Vv cut — -1-1. -1- K-pt^ 


19.74 


0. 


.0696 


2.001 


0.7176 




19.66 


0, 


.0458 


1.575 


5286 


excludes gcil3,xies with lu jC-'^O 


19.72 


0. 


.0625 


1.981 


6326 


iiQpt; nr\ — n 


19.77 


0, 


.0634 


1.965 


0.8923 


uses g3,l3jxies with Si ^0.01 


19.73 


0. 


.0575 


2.135 


0.6276 


iicpQ calaYiPc; 11171+n .S/-">>0 


19.75 


0, 


.0646 


2.108 


0.5981 


uses no selection weighting 


19.73 


0, 


.0637 


2.001 


0.6091 


mean FP shift set to +0.01 


19.73 





.0636 


1.999 


6056 


mpan FP ^Viiff qpI" frt —0 01 


19.71 


0. 


.0639 


1.991 


0.6136 


mpan FP '^Viift '^pt tn -1-0 0'^ 


19.71 


0. 


.0639 


1.992 


0.6104 




19.72 


0, 


.0639 


2.249 


0.4334 


also fit third ajcis of FP 


19.75 


0, 


.0564 


2.192 


0.6402 


uses uniform errors for all galaxies 


19.77 


0, 


.0634 


2.051 


0.6619 


uses clusters with -/Vgg^^ --^3 


19.78 


0, 


.0644 


2.035 


0.6564 


uses clusters with ^^^x --^4 


19.77 


0. 


.0651 


2.084 


0.6139 


iiQPG fliiQtpT't; AS7^il"n /v 

LiOt'O VjL LiO o W 1 U±l 1 V ga,i ^ 


19.74 


0. 


.0643 


1.994 


6060 


hqpc fliicitpT'c; nn^"n J\J i^^f\ 
LiociD ^1 ud o w 1 uii 1 V ga,£ ^ 


19.59 


0, 


.0661 


1.765 


0.5890 


uses clusters with A^gai'^S 


19.52 


0, 


.0544 


1.864 


0.5564 


uses clusters with A^gai^lO 


19.68 


0, 


.0668 


2.009 


0.5470 


excludes spirals 


19.73 


0, 


.0660 


1.966 


0.5488 


excludes spirals; A^gai^3 


19.72 


0, 


.0675 


1.947 


0.5404 


excludes spirals; A^gai^4 


19.72 


0, 


.0678 


2.006 


0.5452 


excludes spirals; A^gai^S 


19.73 


0, 


.0673 


2.027 


0.5414 


excludes spirals; A^gai^G 


19.68 


0, 


.0646 


1.861 


0.6102 


excludes cD galaxies 


19.62 


0, 


.0672 


1.831 


0.5426 


excludes spirals and cDs 



about this FP is (ti=0.064±0.006, corresponding in to an 
intrinsic error in estimating distances of 15%^. Figure |l^ 
sho-ws the projection of the galaxy distribution in the log a- 
FP plane (where FP = r~ bu) . The hard cut in log a and the 
approximate cut in FP are indicated by dashed lines. The 
shape of the best-fitting galaxy distribution is shown by the 
projections of its major and minor axes and its 1, 2, 3 and 4- 
a contours. Figure hp shows the scatter of log Re about the 
FP predictor for log Re , namely alog a + b{SBe) + c. The rms 
scatter about the 1-to-l relation (the solid line) is 0.087 dex, 
which is larger than cri because of the errors in the measure- 
ments. (Allowing for the estimated measurement errors, the 
reduced is 1.01, which is a consistency check on the fit- 
ted value of (Ti.) Thus although the intrinsic rms precision 
of distance estimates from the FP is 0.064 dex (15%), the 
effective rms precision for the EFAR sample when the in- 
trinsic scatter and the measurement errors are combined is 
0.087 dex (20%). 

The random errors given above for the best-fit param- 
eters are based on 1000 simulations of the recovery of the 
FP from the EFAR dataset (assuming no peculiar veloci- 
ties) using the maximum likelihood gaussian algorithm, as 

* Logarithmic errors, e, are converted to linear errors, e, accord- 
ing to e = (10+" - 10-")/2. 



described in Paper VI. Figure g shows the distributions of 
the fitted parameters from these 1000 simulations: the dot- 
ted vertical line is the input value of the parameter and the 
smooth curve is the gaussian with the same mean and rms 
as the fits. There are residual biases in the fitted parameters, 
as shown by the offsets between the input parameters and 
the mean of the fits; a is biased low by 6%, b is biased low 
by 2%, c is biased high by 4%; logi?e, logcrand {S!Be)are 
all biased high, by 0.036 dex, 0.007 dex and 0.05 mag re- 
spectively; the scatter about the FP is under-estimated by 
0.006 dex, or 1.4%. These biases are all less than or compa- 
rable to the rms width of the distribution, so that although 
they are statistically significant (i.e. much greater than the 
standard error in the mean, rms/vTOOO), they do not dom- 
inate the random error in the fitted parameters. We do not 
correct for these biases since they are small and have negli- 
gible impact on the derived distances and peculiar velocities 
(see §H below). 



4.2 Variant cases and systematic errors 

All the other cases listed in Table ^ are variants of this 
standard briefly described in the Notes column of 

Table ^ Case 2 includes galaxies with poorer quality ((5=3) 
photometry and less reliable structural parameters, increas- 
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Table 3. The Fundamental Plane cluster sample. 



CAN 


Name 


cz (kms 


Npp 


X 


A76 


11888 


g 


3 


A119 


13280 


g 


10 


J30 


15546 


g 


13 


A260 


10944 


8 


16 


J8 


9376 


g 


17 


A376 


14355 


7 


20 


A397 


9663 


8 


21 


A400 


7253 


g 


23 


A419 


20400 


g 


24 


A496 


9854 


g 


25 


J34 


11021 


8 


34 


A533 


14488 


g 


35 


A548-1 


11866 


19 


36 


A548-2 


12732 


g 


39 


J13 


8832 


8 


46 


J16W 


11321 


7 


48 


A2040 


13455 


g 


50 


A2063 


10548 


9 


53 


A2147 


10675 


10 


58 


A2151 


11106 


10 


oy 


T1 Q 




7 


65 


A2197 


9137 


9 


66 


A2199 


9014 


9 


68 


A2247 


11547 


7 


70 


J22 


10396 


10 


80 


A2593-N 


12399 


18 


82 


A2634 


9573 


12 


83 


A2657 


12252 


7 


90 


Coma 


6942 


20 



ing the scatter about the FP. Case 3 includes the outlier 
galaxies rejected from the standard sample, and also has a 
larger FP scatter. Cases 4-6 show that applying alternative 
prescriptions for the absorption corrections (BH corrections, 
SFD corrections without an offset, and corrections based on 
a 36:64 weighting of BH and SFD) has no significant effect 
on the FP fit. Case 7 shows that applying a stricter con- 
straint on the errors in the velocity dispersions, excluding 
galaxies for which 51oga>0.1, also has no effect. Cases 8-11 
correspond to different cuts in Dw (no cut and L'vi'cut=6.3, 
14.1 and 15.9 kpc respectively); there is a slight flattening 
of the FP slope a for lower cuts. Case 12 excludes not only 
the galaxies rejected from the standard fit, but also galaxies 
with low likelihoods (ln/I<0); this results in a highly bi- 
ased fit, with both a and b significantly lower than in the 
standard case, and with an artificially lowered FP scatter. 
Case 13 shows that assuming a qo=0 cosmology has no im- 
pact on the FP fit. Cases 14 and 15 examine the effect of a 
lower (5i>0.01) and a higher (Si>0.2) limit on the allowed 
selection probabilities. The former case has highly deviant 
values for logi?e, log a and (SBe) due to over- weighting a 
few galaxies with low selection probabilities; the latter case 
has biased values of a, b and c due to ignoring the tail of the 
selection function. Case 16 ignores the selection probabili- 
ties altogether and applies a uniform weight to all galaxies, 
resulting in an effective over-weighting of the larger galaxies 
and biasing the mean values of log Re and log a to higher 
values. Cases 17-20 show that setting the mean FP shift 
to 4-0.01, -0.01, -1-0.03 and -0.03 dex respectively (rather 
than to zero, as in the standard case) has no effect on the 
fitted FP. 




0.5 1 

a log(c7[km/s]) + b<SB^> + c 



Figure 1. The best-fitting FP solution (case 1) for 255 galaxies 
belonging to the 29 clusters with 6 or more members, (a) The 
projection of the galaxies (marked by their GINs) in the log a- 
FP plane (where FP = r — bu). The dashed lines are the cut in 
logo- and the approximate cut in FP. The best-fitting gaussian 
distribution is shown by the projections of its major and minor 
axes and its 1, 2, 3 and 4-a contours, (b) The scatter of logi?e 
about the FP predictor for log Rf,, namely alog a+b(SBf,)+c. The 
rms scatter about the 1-to-l line is 0.087 dex (an rms distance 
error of 20% per galaxy). The inset histogram of residuals Alog Rg 
has a gaussian with an rms of 0.087 dex overlaid. 



Case 21 permits an extra degree of freedom by allowing 
the orientation of the major axis of the galaxy distribution 
within the FP to be fitted, rather than specified a priori. 
The unit vectors of the galaxy distribution for the standard 
case, given by equation (H) , are 
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1.223 -• 1.150±0.0 



: 0-336 -> 329±0.013 



: -8.666 -• -8.360±0.351 




0.25 0.3 0.35 




1.1 1.2 1.3 1.4 



-7.5 








■12 




-8 






-8.5 


- " ^'^ ^ 




-9 







0.3 0.32 0.34 0.36 



<lagR,> ; 0.770 0.808±0.025 



log(7> : 2.304 -. 2.311±0.016 <SB,> : 19.71 -> 19.77±0.06 




0.04 0.06 0. 



Figure 2. The distributions of the FP parameters a, b, c, logiJe, 
logcr, (SBe)i Ti, cr2 and (T3 resulting from fitting 1000 simulations 
of the best-fit FP. The input parameters of the simulations arc 
given at the head of each panel (and indicated by the vertical 
dotted line), followed by the mean and rms of the fits to the 
simulations (the curve is the gaussian with this mean and rms). 



vi = +1.000r - 1.223s - 0.336m 
V2 = +1.000f + 0.000.s + 2.978m 
V3 = -2.978f - 2.710s + l.OOOii 



(11) 



while the true eigenvectors, obtained by fitting with the ex- 
tra degree of freedom, are 



VI = -l-l.OOOf - 1.227s - 0.337U 
V2 = +1.000f - 0.032s + 2.964U 
V3 = -3. 176f - 2.863s + l.OOOu 



(12) 



The coefficient of s in the second eigenvector is small, justi- 
fying the simplifying approximation of setting it to zero used 
in equation (W). The FP itself is very close to the standard 
fit, while the axes within the FP have coefficients differing 
from the standard values by no more than a few percent; ai 
stays the same, while (T2 is maximised and 0-3 is minimised. 

Case 22 replaces the individual error estimates for all 
measured quantities with uniform errors; this has little effect 
on the FP, but under-estimates the intrinsic scatter about 
the plane. Cases 23-28 explore the effects of varying the 
minimum number of galaxies required for a cluster to be in- 
cluded in the fit, from 3, 4, 5, 6 and 8 up to 10. Note that 
this is the number of galaxies in the cluster before excluding 
outliers; hence case 26 differs from case 1 in having 31 clus- 
ters rather than 29. The simulations of Paper VI suggested 
that a spuriously small estimate for cti could in principle 
result when clusters with few galaxies are included in the 
fit, as offsetting the FP with a spurious peculiar velocity 
could suppress the apparent scatter. However this effect is 
not observed in fitting the actual data, and the FP fits are 
consistent with the errors on the best fit for samples with a 
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Figure 3. The fitted FP parameters for each case in Table ^ 
showing the distributions and correlations for various pairs of pa- 
rameters. Each case is numbered as in the table. The dots are the 
distribution of fits obtained for 1000 simulations of the standard 
case after removing the effects of the residual biases. 



minimum number of galaxies per cluster of between 3 and 
8. A significantly fiatter FP slope is found only for the set 
of clusters with 10 or more galaxies, where there are only 7 
clusters and 99 galaxies in the fit and correspondingly larger 
uncertainties. Case 29 is the same as the standard case ex- 
cept that spirals are excluded, so that the FP is fitted only to 
galaxies with E, E/SO and cD morphological types. The FP 
slope for these early-type galaxies is steeper, with a=1.33. 
Cases 30-33 are similar to case 29, but with the minimum 
number of galaxies required for a cluster to be included in 
the fit varied from 3 to 6. Cases 34 and 35 are the same as the 
standard case except that the fit is restricted, respectively, 
to exclude cD galaxies and both cD galaxies and spirals. Re- 
moving cDs flattens a and lowers b, in contrast to case 29; 
removing both cDs and spirals restores the FP to the inter- 
mediate values obtained by including both populations. 

Figure ^ shows the fitted values in each case for various 
pairs of the parameters, in order to show their distributions 
and correlations. The cases are numbered following Table ^ 
The dots show the distribution of fits obtained for 1000 sim- 
ulations of the standard case (case 1) after removing the ef- 
fects of the residual biases. The main point to note is that, 
with only a few exceptions (noted above), the systematic 
differences in the fits derived for difference cases are compa- 
rable to the random errors in the determination of the pa- 
rameters for the standard case. We conclude that the uncer- 
tainties in our best-fit FP parameters are dominated by the 
random errors and not by systematic effects from the fitting 
method. In particular we conclude that the following inputs 
have relatively little effect on the fitted FP: the absorption 
correction, the cosmological model, the assumed mean FP 
shift and the choice of the second and third FP axes. Our 
standard case provides an optimum fit to the FP because: 
(i) it excludes the galaxies with poor structural parameters 
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Table 5. Determinations of the Fundamental Plane. 



Source 


Band 


A 


B 


A 


Fit method 


Dressier et al. (1987) 


B 


1.33 ± 0.05 


—0.83 


i 0.03 


20% 


inverse 


Djorgovski & Davis (1987) 




1.39 ± 0.14 




-1- n riQ 


20% 


2-step inverse 


Lucey et al. (1991) 


B 


1.27 ± 0.07 


-0.78 


± 0.09 


13% 


inverse 


Guzman et al. (1993) 


V 


1.14 


-0.79 




17% 


forward 


j0rgensen et al. (1996) 


r 


1.24 ± 0.07 


-0.82 


± 0.02 


17% 


orthogonal 


Hudson et al. (1997) 


R 


1.38 ± 0.04 


-0.82 


± 0.03 


20% 


inverse 


Scodeggio et al. (1997) 


I 


1.25 ± 0.02 


-0.79 


± 0.03 


20% 


mean regression 


Pahre et al. (1998) 


K 


1.53 ± 0.08 


-0.79 


± 0.03 


21% 


orthogonal 


Miiller et al. (1998) 


R 


1.25 


-0.87 




19% 


orthogonal 


Gibbons et al. (2000) 


R 


1.39 ± 0.04 


-0.84 


± 0.01 


19% 


inverse 


EFAR (this paper) 


R 


1.22 ± 0.09 


-0.84 


± 0.03 


20% 


ML gaussian 



and velocity dispersion measurements which artificially in- 
flate the scatter about the FP and the uncertainty in the 
FP parameters; (ii) it applies a selection function cutofi' that 
balances over-weighting a small number of galaxies against 
biasing the results by ignoring galaxies with low selection 
probabilities; (iii) it uses clusters with 6 or more galaxies 
to avoid artificially reducing FP scatter by confusing scat- 
ter with peculiar velocities while yet retaining a sufficiently 
large overall number of galaxies to keep the random errors 
in the FP parameters small. 



4.3 Comparison with previous work 

Table ^ compares the best-fit EFAR FP with previous de- 
terminations in the literature, noting both the passband to 
which the relation applies and the method of the fit. To 
match the usage in most of this literature, we present the 
FP in the form Re oc (t^{E)^, where ctq is the central ve- 
locity dispersion and {E)e is the mean surface brightness (in 
linear units) within the effective radius Re- The exponents 
of this relation are related to the coefficients of our FP re- 
lation, log Re = alog a + b(SBe) + c, by A— a and B=—2.5b. 
The table also quotes the fractional distance error. A, cor- 
responding to the rms scatter about the FP in Re- In most 
cases the determination of the FP is limited to galaxies with 
cr>100kms~^. The forward and inverse fitting methods are 
linear regressions with, respectively, log Re and log a as the 
independent variable; orthogonal fitting minimises the resid- 
uals orthogonal to the FP, while mean regression averages 
the fits obtained by taking each of logi?e, logcr and (SBe) 
as the independent variable. 

The first point to note is that all the fitted values of B 
are consistent within the errors, regardless of passband and 
fitting method. The second point to note is that this is not 
true for A, which has a higher value in the A'-band FP fit 
of Pahre et al. (1998) than in any of the optical fits. The 
third point is that, within the optical FP fits, the forward 
and inverse fits give, respectively, lower and higher values of 
A than the orthogonal and mean regressions and the maxi- 
mum likelihood gaussian method. This is consistent with the 
analysis and simulations of the methods carried out in Pa- 
per VI: for samples in which the errors in a dominate and/or 
selection cuts are applied in Re (as is the case for most of 
these datasets), the value of A will be under-estimated by a 
forward fit and over-estimated by an inverse fit. Orthogonal 
and mean regressions reduce these biases, with the least bias 
being produced with the maximum likelihood method. We 



conclude that apparent differences between FP fits in opti- 
cal passbands are due to differences in the fitting methods 
that have been applied. 

There is also consistency on the observed scatter about 
the FP as represented by the fractional distance error, A. 
With the exception of Lucey et al. (1991), the observed er- 
rors are all in the range 17% to 21%. This is consistent with 
(i.e. larger than) the estimated intrinsic scatter about the 
FP of 15% that we derive from the EFAR sample, and the 
range corresponds to the range of measurement errors in the 
various studies, which account for between 8% and 15% of 
the observed scatter. 



5 DISTANCES AND PECULIAR VELOCITIES 

In order to determine distances and peculiar velocities, we 
re-apply the maximum likelihood gaussian algorithm to the 
whole cluster sample. This time we fix the parameters of the 
intrinsic galaxy distribution at their best-fit values (case 1 of 
Table m and fit only for the shift of the FP for each cluster. 



5.1 Sample 

We remove outliers (interlopers from the cluster foreground 
or background, objects which genuinely do not lie on the FP, 
and objects with bad data) by excluding the galaxies that 
deviate most from the fitted FP until all clusters have FP fits 
with x^/i'<3. To check that this procedure is conservative, 
we visually inspected each cluster's distribution of Dg — De 
(individual galaxy distances relative to the overall cluster 
distance, from the residuals about the best-fit FP) with re- 
spect to czg — czc (individual galaxy redshifts relative to the 
overall cluster redshift). The rejected galaxies were invari- 
ably clear outliers in these distributions. In all, 36 galaxies 
were rejected using this procedure, including all the galaxies 
rejected from the FP fit (see Table ^). The list of galaxies 
excluded from the peculiar velocity fits is given in Table ^ 
There were three clusters with x^/z^>3 (CAN 2=A85 with 4 
galaxies, CAN 55=P386-2 with 2 galaxies, CAN 79=A2589 
with 5 galaxies) for which half or more of the galaxies had 
to be rejected in order to obtain a good FP fit, so that it 
was difficult to determine which galaxies were the outliers. 
Although we give distances and peculiar velocities for these 
clusters below (using all the available galaxies), we exclude 
them from further analysis. 
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Table 6. Galaxies excluded from the peculiar velocity fits. 
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25 
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66 
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25 


J34 F 
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80 


A2593-S C 
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35 


A548-1 F 


713 


80 


A2593-S E 


355 


42 


J14-1 D 


721 


82 


A2634 F 


370 


43 


A1983 2 


728 


82 


A2634 2 


396 


46 


J16-W B 


730 


83 


A2657 B 


432 


50 


A2063 G 


731 


83 


A2657 C 


456 


53 


A2147 D 


756 


90 


COMA 133 



5.2 Bias corrections 

To the extent that its assumptions are justified, the maxi- 
mum Ukelihood gaussian algorithm accounts for the efi'ects 
of biases on the estimated distances which are due to the se- 
lection function of the galaxies within each cluster. (We refer 
to this bias as 'selection bias' rather than 'Malmquist bias' 
because, following the usage of Strauss & Willick (1995), the 
effect is due to the selection criteria and not the line-of-sight 
density distribution.) As discussed in Paper VI, however, the 
sample selection function parameter varies with cluster 
redshift, introducing a redshift-dependent bias in the pecu- 
liar velocity estimates. Although this bias is reduced by the 
selection limit Dwi>Dwcut imposed on galaxy sizes (see 

clusters with D'^j>Dwcut are nonetheless sampled dif- 
ferently to clusters with D^Yj<Dwcut- This difference in the 
way the FP galaxy distribution is sampled in different clus- 
ters leads to a residual bias in the clusters' fitted FP offsets 
and peculiar velocities as a function of D^rj (or redshift, 
with which is closely correlated). 

This effect is investigated in detail through simulations 
in Paper VI. Figure ^ shows the residual selection bias de- 
termined from 1000 simulations of the EFAR dataset. For 
clusters with redshifts below the sample mean the bias in 
the peculiar velocities is small and negative, while for clus- 
ters at redshifts above the sample mean it is positive and 
increases rapidly with redshift. We correct this systematic 
bias individually for each cluster by subtracting the mean 
error in the FP offset determined from 1000 simulations of 
the EFAR dataset before computing the cluster distances 
and peculiar velocities. The size of the corrections are shown 
in the inset histograms of Figure ^. For the subsample of 
clusters included in subsequent analyses of the peculiar ve- 
locities (whose selection is discussed below), the amplitude 
of the bias correction is less than 250kms~^ for 40 of the 
50 clusters. The random errors in the peculiar velocities are 
typically of order 1000 km s~^, while the uncertainties in the 
peculiar velocity bias corrections for these clusters are typ- 
ically less than 50kms~^. To the extent that the simulated 
datasets match the real distribution of galaxies in the FP, 
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Figure 4. The residual selection bias determined from 1000 sim- 
ulations of the EFAR dataset. (a) The bias in the FP offsets {Sj) 
for each cluster as a function of the cluster's selection function 
parameter f^j- The inset histogram shows the distribution of 
bias corrections (S) . (b) The corresponding bias in the cluster pe- 
culiar velocities (Vpcc) as a function of cluster redshift cz. The 
inset histogram shows the distribution of bias corrections (Vpcc). 
The filled symbols and the shaded histogram show the subsample 
of clusters used in the peculiar velocity analysis. 



therefore, the bias corrections should not significantly in- 
crease the random errors in the peculiar velocities. 



5.3 Results 

The individual FP fits are shown in Figure where the fixed 
parameters of the galaxy distribution used to fit the FP shift 
are given at the top of the plot. Each panel corresponds to 
a cluster, labelled by its CAN; the 29 clusters used to derive 
the parameters of the galaxy distribution are indicated by 
bold labels. The area of each point is proportional to the 
selection weight of the galaxy; the corresponding GINs are 
given at left. The solid line is the major axis of the global fit 
to the FP, and the cross on this line the centre of the global 
galaxy distribution, (log a, log Re—b(SBe)). The dotted lines 
and ellipse are the major and minor axes and the 3a contour 
of the cluster's FP, vertically offset from the global FP by the 
cluster's FP shift. The cluster's mean redshift cz, distance 
D, and peculiar velocity Vpec, each with its estimated error, 
are given at the bottom of the panel. The distances and 
peculiar velocities are corrected for the residual selection 
bias discussed above. 
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Figure 5. The FP distributions for each individual cluster. See text for description. 
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Figure 5. Continued. 
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Figure 5. Continued. 
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Figure 5. Continued. 
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Figure 5. Continued. 



The results are summarised in Table |7| which for each 
cluster gives: CAN, cluster name (in parentheses for fore- 
and background groups), number of galaxies used in the 
distance determination. Galactic longitude and latitude, the 
bias-corrected FP shift 5 and its uncertainty, the bias cor- 
rection es that was subtracted from the raw value of 5, the 
cluster redshift cz and its uncertainty Acz, and the bias- 
corrected values of the cluster distance D and its uncertainty 
AD, the redshift czd corresponding to D, and the peculiar 
velocity V and its uncertainty AV . Note that some clusters 
are missing from this list; CAN 81 because it has been com- 
bined with CAN 80 (see Paper II); CANs 41, 47, 54 and a 
number of the fore- and background groups (CANs > 100) 
because no cluster members meet the selection criteria. 



5.4 Tests and comparisons 

Gibbons et al. (2000) have suggested that the large peculiar 
velocities found for some clusters are due to poor FP fits. 
For a heterogeneous sample of 20 clusters drawn from their 
own observations and the literature, they find that nearly 



half are poorly-fit by a FP and have twice the rms scatter of 
the well-fit clusters. The half of their clusters that have good 
FP fits all have peculiar velocities that are consistent with 
them being at rest in the CMB frame; the poorly-fit clus- 
ters show a much larger range of peculiar velocities. Gibbons 
et al. suggest that the large peculiar velocities detected for 
some clusters may result from those clusters being poorly 
fit (for whatever reason) by the global FP. The origin of the 
poor fits is not known, but the possibilities include intrinsic 
FP variations between clusters, failure to identify and re- 
move interlopers, observational errors, the heterogeneity of 
the data, and combinations of these efi'ects. 

We therefore need to test whether some of the pecu- 
liar velocities we derive from the EFAR dataset are due to 
poor fits to the FP rather than genuine peculiar velocities. 
Figure ^ shows the peculiar velocities of the EFAR clus- 
ters as a function of the goodness-of-fit of their best-fit FP 
(as measured by the reduced statistic). As noted above, 
even after removing outliers, there are still three clusters 
with very poor FP fits (x^/j^>3; in fact CAN 2=A85 actu- 
ally has /v—11, but is plotted at /v=^ for convenience). 
All three of these clusters have large negative peculiar veloc- 
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Peculiar motions of early-type galaxies - VII 

Table 7. Redshifts, distances and peculiar velocities for the EFAR clusters in the CMB frame 
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+47.23 


-0.0736 


0.0887 


-0.0316 


26359 


339 


19195 


3620 


21696 


4643 


+4348 


4364 


57* 


J18 


3 


39.95 


+46.50 


-0.0743 


0.0591 


-0.0032 


9661 


132 


7702 


1045 


8069 


1147 


+1549 


1124 


58* 


A2151 


10 


31.47 


+44.64 


-0.0469 


0.0520 


-0.0165 


11194 


79 


9422 


1085 


9978 


1218 


+1176 


1182 


59* 


J19 


5 


29.06 


+43.50 


-0.0966 


0.0351 


-0.0047 


12782 


215 


9510 


761 


10077 


855 


+2616 


854 


60* 


P445-1 


5 


31.19 


+46.17 


-0.0057 


0.0374 


-0.0251 


13804 


119 


12591 


1015 


13607 


1188 


+187 


1145 


61 


P445-2 


2 


28.77 


+45.63 


-0.0632 


0.0499 


+0.0083 


4978 


90 


4181 


491 


4288 


517 


+680 


517 
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Table 7. Continued. 



CAN 


Name 




I 


b 


5 


A(5 




cz 


Acz 


D 


AD 


C2£3 


Aczo 


V 


AV 


62 


A2162-N 


5 


50.36 


+46.10 


-0.0629 


0.0499 


-0.0066 


15048 


154 


11948 


1352 


12859 


1569 


+2098 


1513 


63 


A2162-S 


1 


48.36 


+46.03 


-0.0892 


0.0872 


+0.0010 


9837 


116 


7570 


1537 


7925 


1685 


+ 1862 


1646 


64 


J20 


2 


56.54 


+45.58 


+0.0221 


0.0710 


-0.0309 


9468 


113 


9434 


1433 


9991 


1608 


-506 


1564 


65* 


A2197 


7 


64.68 


+43.50 


-0.0566 


0.0435 


—0.0015 


9159 


97 


7627 


762 


7987 


836 


+1141 


820 


66* 


A2199 


7 


62.92 


+43.70 


+0.0198 


0.0320 


+0.0001 


9039 


93 


8980 


664 


9483 


741 


—431 


723 


67* 


J21 


4 


77.51 


+41.64 


+0.0174 


0.0419 


-0.0198 


13876 


215 


13343 


1224 


14492 


1447 


-587 


1398 


68* 


A2247 


7 


114.45 


+31.01 


-0.0127 


0.0576 


-0.0060 


11506 


67 


10460 


1371 


11150 


1560 


+342 


1506 


69 


P332-1 


1 


49.95 


+35.22 


-0.1640 


0.0874 


—0.0217 


17338 


403 


10770 


2067 


11503 


2361 


+5619 


2311 


70* 


J22 


10 


49.02 


+35.93 


-0.0581 


0.0319 


—0.0026 


10396 


78 


8566 


625 


9023 


694 


+1332 


678 


71* 


J23 


4 


85.81 


+35.40 


+0.0535 


0.0501 


+0.0025 


8401 


113 


9052 


1054 


9564 


1177 


— 1127 


1146 


72* 


J24 


3 


69.51 


+32.08 


-0.0088 


0.0722 


—0.0096 


10256 


136 


9476 


1547 


10038 


1737 


+210 


1687 


73 


J25 


2 


91.82 


+30.22 


-0.0565 


0.0661 


+0.0013 


7863 


194 


6598 


1012 


6865 


1096 


+974 


1088 


74 


J26 


2 


69.59 


+26.60 


-0.0728 


0.0680 


-0.0329 


14876 


267 


11556 


1668 


12405 


1925 


+2372 


1872 


75* 


J27 


4 


80.41 


+23.15 


-0.0084 


0.0541 


+0.0001 


9667 


233 


8968 


1121 


9470 


1251 


+190 


1234 


76 


J38 


1 


36.09 


-44.90 


-0.0256 


0.0887 


—0.0306 


15056 


320 


13027 


2476 


14119 


2914 


+894 


2809 


77* 


P522-1 


3 


81.75 


-41.26 


-0.0770 


0.0393 


+0.0026 


7362 


112 


5910 


539 


6124 


579 


+ 1212 


578 


78* 


A2572 


4 


94.28 


-38.95 


+0.0444 


0.0374 


—0.0075 


11210 


136 


11643 


984 


12505 


1137 


— 1243 


1100 


79 


A2589 


5 


94.64 


-41.23 


+0.0751 


0.0309 


-0.0022 


12078 


149 


13399 


947 


14558 


1121 


-2365 


1079 


80* 


A2593-N 


16 


93.44 


-43.19 


+0.0328 


0.0170 


—0.0016 


12034 


107 


12113 


472 


13051 


549 


-974 


537 


82* 


A2634 


10 


103.50 


-33.08 


-0.0026 


0.0225 


-0.0021 


9228 


113 


8698 


448 


9169 


498 


+56 


495 


83* 


A2657 


5 


96.73 


-50.25 


+0.0644 


0.0396 


—0.0024 


11887 


156 


12879 


1170 


13945 


1375 


— 1966 


1322 


84 


A2666 


2 


106.71 


-33.80 


-0.0365 


0.0545 


+0.0035 


7620 


90 


6704 


852 


6981 


924 


+624 


907 


90* 


COMA 


19 


58.00 


+88.00 


+0.0017 


0.0173 


+0.0018 


7209 


50 


6942 


277 


7238 


302 


-29 


299 


108 


(A168) 


2 


134.36 


-61.61 


+0.0861 


0.0532 


+0.0055 


4683 


136 


5556 


692 


5745 


740 


-1042 


738 


109 


(A189) 


1 


140.13 


-59.99 


-0.0118 


0.0851 


—0.0191 


9165 


126 


8461 


1592 


8906 


1765 


+251 


1721 


120 


(A397) 


1 


161.84 


-37.33 


+0.0402 


0.0709 


-0.0314 


16631 


15 


16596 


2499 


18423 


3089 


— 1688 


2923 


128 


(J35) 


1 


217.47 


-33.61 


-0.0359 


0.0812 


+0.0000 


16907 


27 


14140 


2668 


15439 


3188 


+ 1395 


3032 


129 


(J34/35) 


1 


216.40 


-34.19 


-0.1005 


0.0830 


-0.0466 


16464 


18 


11896 


2027 


12798 


2350 


+3515 


2264 


130 


(P777-1) 


1 


218.49 


-32.70 


+0.0377 


0.0746 


-0.0240 


16755 


14 


16615 


2686 


18447 


3322 


-1594 


3141 


131 


(P777-2) 


1 


220.77 


-32.62 


-0.0834 


0.0698 


+0.0007 


9300 


102 


7275 


1178 


7602 


1286 


+ 1655 


1258 


132 


(P777-3) 


2 


219.72 


-31.71 


-0.0070 


0.0482 


+0.0000 


13502 


144 


12300 


1369 


13268 


1596 


+224 


1534 


136 


(A548-2) 


1 


230.40 


-25.97 


-0.0031 


0.0706 


-0.0011 


8826 


234 


8329 


1358 


8760 


1503 


+63 


1478 


138 


(J12) 


1 


50.52 


+78.23 


-0.0580 


0.0740 


-0.0124 


15633 


16 


12512 


2077 


13515 


2428 


+2026 


2327 


142 


(J14-1) 


1 


8.80 


+58.73 


-0.0625 


0.0936 


-0.0490 


16121 


125 


12736 


2405 


13778 


2819 


+2239 


2711 


144 


(A1991) 


1 


22.74 


+60.52 


+0.1807 


0.1143 


-0.0374 


13372 


23 


18781 


4529 


21166 


5774 


—7280 


5426 


145 




3 


6.81 


+48.20 


-0.0985 


0.0553 


-0.0119 


15927 


109 


11594 


1437 


12449 


1659 


+3338 


1598 


146 


(J16W) 


1 


5.08 


+49.63 


-0.1752 


0.0870 


+0.0000 


24347 


221 


14186 


2871 


15494 


3433 


+8417 


3271 


154 


(P386-1) 


1 


37.09 


+47.81 


-0.0679 


0.1159 


—0.0129 


13910 


212 


10988 


2889 


11752 


3309 


+2075 


3195 


160 


(P445-1) 


1 


31.19 


+46.17 


+0.1585 


0.0783 


-0.0134 


10827 


223 


14657 


2569 


16059 


3092 


-4966 


2948 


163 


(A2162-S) 


3 


48.36 


+46.03 


+0.0706 


0.0750 


-0.0138 


15623 


273 


16817 


2817 


18697 


3493 


—2894 


3305 


166 


(A2199) 


1 


62.92 


+43.70 


-0.0371 


0.1774 


-0.0719 


17820 


44 


14785 


5183 


16213 


6248 


+1524 


5978 


167 


(J21) 


1 


77.51 


+41.64 


-0.0867 


0.0852 


+0.0025 


6086 


12 


4812 


956 


4953 


1013 


+ 1113 


997 


174 


(J26) 


1 


69.59 


+26.60 


+0.0896 


0.1401 


-0.0788 


18098 


36 


20073 


5312 


22828 


6899 


-4396 


6496 


177 


(P522-1) 


1 


81.75 


-41.26 


+0.0033 


0.0710 


-0.0300 


23626 


66 


20842 


3145 


23832 


4130 


-191 


3848 


180 


(A2593-N) 


1 


93.44 


-43.19 


-0.1306 


0.0696 


-0.0324 


27458 


8 


17433 


2570 


19464 


3215 


+7505 


3034 


210 


(J30) 


1 


151.84 


-75.04 


+0.0180 


0.0676 


—0.0234 


20667 


12 


19167 


2805 


21659 


3596 


-925 


3368 


232 


(P777-3) 


1 


219.72 


-31.71 


+0.0855 


0.0673 


+0.0009 


8929 


9 


10327 


1612 


10999 


1831 


-1997 


1766 


238 




1 


50.52 


+78.23 


-0.1249 


0.1134 


-0.0422 


18073 


27 


12235 


2904 


13192 


3382 


+4674 


3253 


239 


(J13) 


1 


28.27 


+75.54 


-0.0540 


0.1141 


-0.0477 


17675 


5 


14118 


3316 


15412 


3961 


+2151 


3788 


240 


(J36) 


1 


332.77 


+49.31 


+0.0445 


0.0816 


-0.0284 


17559 


26 


17607 


3081 


19681 


3863 


-1991 


3641 


244 


(A1991) 


1 


22.74 


+60.52 


+0.0315 


0.1564 


-0.1206 


21381 


64 


20373 


5377 


23219 


7014 


-1706 


6638 


245 


(J16) 


1 


6.81 


+48.20 


-0.1560 


0.1049 


-0.0198 


14194 


39 


9138 


2126 


9660 


2377 


+4392 


2306 


254 


(P386-1) 


1 


37.09 


+47.81 


+0.0711 


0.2459 


-0.1072 


26704 


73 


27084 


12329 


32415 


17745 


-5154 


16417 


260 


(P445-1) 


1 


31.19 


+46.17 


+0.1227 


0.0970 


-0.0630 


18642 


26 


22247 


4181 


25696 


5604 


-6497 


5225 


338 


(J12) 


1 


50.52 


+78.23 


-0.0493 


0.0933 


-0.0066 


10666 


236 


8955 


1911 


9455 


2132 


+1173 


2081 


339 


(J13) 


1 


28.27 


+75.54 


+0.0993 


0.1160 


-0.0367 


15883 


30 


18235 


4476 


20472 


5662 


-4296 


5331 


340 


(J36) 


1 


332.77 


+49.31 


-0.1612 


0.1989 


+0.0000 


29241 


28 


17136 


8373 


19093 


10428 


+9540 


9804 



Note: clusters in the Fundamental Plane sample (Table ^) have their CANs in bold; clusters in the peculiar velocity sample are 
marked with an asterisk. 



This table is also available as J /MNKAS/ vol /page from NASA's Astrophvsical Data Centre (ADC. |ittp:/ /adc. gsfc.nasa.gov ) and 
from the Centre de Donnees astronomiques de Strasbourg (CDS, tittp;//cdsweb.u-strasbg.fr ). 
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Figure 6. Cluster peculiar velocities as a function of the 
goodness-of-fit of their best-fit FP. Clusters indicated by their 
CANs; those with only single members have no x^/^ ^i^d are 
plotted at the left of the figure to show their peculiar velocities. 
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Figure 7. The distribution of cluster peculiar velocities with re- 
spect to the cluster offsets from the global Mg— cr relation of Pa- 
per V. 



ities, detected at nominal significance levels of 1.8-2.5(t. The 
poor quality of the FP fits raises considerable doubts about 
the reality of the peculiar velocity estimates, however, and 
we therefore omit these clusters from all subsequent anal- 
ysis. The remaining clusters generally have acceptable fits 
(x^/u^l). There are 10 clusters with x^/i'=2-3, but none 
of these have significant peculiar velocities (the strongest 
detection is at the l.Scr level). Apart from the three clusters 
with x^/j'>3, the clusters are all adequately fitted by the 
global FP, and there is no evidence for any increased scatter 
in the peculiar velocities for poorer FP fits. 

Another possible source of systematic errors are the 
small biases in the recovered parameters of the best-fitting 
FP (see ^ above). If we apply the corrections for these bi- 
ases derived from our simulations (Figure ^ and re-derive 
the peculiar velocities with this bias-corrected FP, we find 
that the peculiar velocities of the clusters are not signif- 
icantly altered: the peculiar velocity of Coma changes by 
-fl4kms~^, and the rms difference in peculiar velocity be- 
tween our standard solution and the bias-corrected solution 
is only 67kms~^. 

We can also attempt to test whether differences in the 
mean stellar populations between clusters produce spurious 
peculiar velocities, by looking for a correlation between the 
peculiar velocities and the offset of each cluster from the 
global Mg-cr relation derived in Paper V. The correlation 
coefficient for the distribution (shown in Figure ^ is —0.30, 
but 1000 simulations of the observed distribution show that, 
allowing for the estimated errors, this value does not indicate 
a correlation significant at the 95% level. However, while 
there is no positive evidence that stellar population differ- 
ences are leading to spurious peculiar velocities, this test 
cannot rule out this possibility. Figure 10 of Paper V shows 
that the joint distribution of residuals about the FP and 
Mg-(j relations is consistent with simple stellar population 



models if one invokes sufficiently large (and possibly corre- 
lated) scatter in the ages and metallicities of the galaxies. 
Against this possibility we can set the generally good agree- 
ment between the distance estimates obtained from the FP 
and other methods (such as the TuUy-Fisher relation and 
surface brightness fluctuations) which have different depen- 
dences on the stellar populations. 

Finally, we can perform a direct comparison between 
the peculiar velocities we measure and those obtained by 
other groups for the same clusters. Figure ^ shows com- 
parisons with the TuUy-Fisher estimates of Giovanelli et al. 
(1998b; SCI) and Dale et al. (1999b; SCII), and the FP es- 
timates of Hudson et al. (1997; SMAC) and Gibbons et al. 
(2000; GFB). The flattening in the VeFAR-Vothor distribu- 
tions is due to the fact that the uncertainties in the EFAR 
peculiar velocities are generally larger than those of the 
other measurements — although the error per galaxy is sim- 
ilar in all cases, the EFAR sample typically has a smaller 
number of galaxies per cluster. A x^-test shows that the pe- 
culiar velocity measurements are consistent within the errors 
in all three comparisons. 



6 BULK MOTIONS 
6.1 Cluster sample 

In analysing the peculiar motions of the clusters in the 
EFAR sample we confine ourselves to the subsample of 
clusters with 3 or more galaxies {Ng>3), C2:< 15000 km s~^ 
and (Sy<1800kms~^. These criteria are illustrated in Fig- 
ure 1^, and are chosen because: (i) they eliminate all the 
fore- and background clusters, for which the selection func- 
tions have not been directly measured and are only poorly 
approximated by the selection function of the main clus- 
ter onto which they are projected; (ii) they eliminate the 
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Figure 8. Comparisons of EFAR peculiar velocities of clusters in 
common with SCI/II (Giovanelli et al. 1998b, Dale et al. 1999b), 
SMAC (Hudson et al. 1997) and GFB (Gibbons et al. 2000). 



clusters with only 1 or 2 galaxies in the FP fit, where it is 
not possible to check if galaxies are cluster interlopers or 
FP outliers; (iii) they eliminate the higher-redshift clusters, 
which have proportionally higher uncertainties in their pe- 
culiar velocities (and in any case sample the volume beyond 
c«= 15000 km s~^ too sparsely to be useful); (iv) they elimi- 
nate clusters with large uncertainties in their peculiar veloc- 
ities, resulting from large measurement errors for individual 
galaxies exacerbated by a small number of galaxies in the 
cluster — restricting the subsample to (5V<1800kms~^ (the 
peculiar velocity error for a cluster with a FP distance from 3 
galaxies with a distance error per galaxy of 20%) represents a 
compromise between using clusters with better-determined 
peculiar velocities and keeping the largest possible cluster 
sample. 

We also eliminate from the sample the three clusters 
which were identified in the previous section as having un- 
acceptably poor FP fits (CAN 2=A85, CAN 55=P386-2 and 
CAN 79=A2589); two of these would be eliminated in any 
case: A85 because it has C2>15000 kms~^, and P386-2 be- 
cause it has only two galaxies. We also eliminate the two 
components of A548 (CAN 35=A548-1 and CAN 36=A548- 
2), since the substructure in this region (ZabludofT et al. 
1993, Davis et al. 1995) makes cluster membership prob- 
lematic and since the high relative velocity of the two main 
subclusters is not relevant to the large-scale motions we are 
investigating (Watkins 1997). 

The subsample selected in this way for the analysis of 
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5000 10000 15000 

cz (km/s) 

Figure 9. The selection of clusters for the peculiar velocity analy- 
sis. The cluster's peculiar velocity errors are plotted as a function 
of their redshifts. Each cluster is marked by its CAN, with clusters 
having 3 or more galaxies in a larger font; fore- and background 
groups (CAN>100) are not shown. The selection limits in cz and 
SV are indicated by the dotted lines. The distribution of peculiar 
velocity errors is shown in the inset: the open histogram is for all 
clusters, the filled histogram for the selected clusters. 



the peculiar motions comprises 50 clusters (25 in HCB, 25 
in PPC); they are indicated by an asterisk in Table ^. The 
distribution of the peculiar velocity uncertainties for this 
subsample is shown in the inset to Figure |^); the median 
peculiar velocity error is 1060kms"\ Figure ^ shows the 
projection of the sample on the sky in Galactic coordinates, 
with the amplitude of the clusters' peculiar velocities in the 
CMB frame indicated by the size of the symbols. Inflow- 
ing clusters (circles) and outflowing clusters (asterisks) are 
fairly evenly distributed over the survey regions. The median 
direction of the clusters belonging to the peculiar velocity 
sample in the HCB region is (/,&)=(42°,48°), and in the PPC 
region is (Z,6) = (152°,-36°); the angle between these two di- 
rections is 128°. 



6.2 Bulk motions 

The peculiar velocities of the sample clusters as a function 
of redshift are shown in Figure [111. The mean peculiar ve- 
locity of the whole sample ((V) = 159±158kms~^) is consis- 
tent, within the errors, with no net inflow or outflow. This 
need not have been the case, as the FP zeropoint is based 
on the 29 clusters listed in Table |^, which make up only 
26 of the 50 clusters in the peculiar velocity sample. The 
mean peculiar velocities of each of the two sample regions 
separately are also consistent with zero inflow or outflow: 
(l^HCB)=+383±229kms-^ {l^ppc)=-65±217kms-\ A 
test shows that the observed peculiar velocities are consis- 
tent with strictly zero motions (i.e. no bulk or random mo- 
tions at all) at the 2% level. If the one cluster with a 3(j pecu- 
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Figure 10. The projection on the sky in Galactic coordinates of 
the EFAR peculiar velocities in the CMB frame. Clusters with 
positive (negative) peculiar velocities are indicated by asterisks 
(circles) ; marker sizes are related to the amplitude of the peculiar 
velocity. Other markers show the directions with respect to the 
CMB frame of the Local Group dipole (0), the Lauer &; Postman 
(1994) dipole (0), the SMAC (Hudson et al. 1999) dipole (0), 
and the LPIOK (Willick 1999) dipole (0). 
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Figure 12. The mean peculiar velocity in radial shells. The 
clusters are grouped into 7 redshift ranges: the first is 4000- 
8000kms-l, the next five cover 8000 kms"! to 13000 kms"! in 
lOOOkms-l steps, and the last is 13000-15000 km s"!. The left 
panel shows the whole sample of 50 clusters, the middle panel 
shows the 25 HCB clusters, and the right panel shows the 25 
PPC clusters. The Supergalactic X, Y and Z components are 
shown as filled squares, circles and triangles respectively (with 
small offsets in redshift for clarity). The number of clusters in 
each redshift range is indicated at the bottom of each panel. 
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Figure 11. Peculiar velocities of the EFAR clusters as a function 
of redshift. The clusters in PPC are given negative redshifts, but 
in all cases positive peculiar velocities indicate outflow and nega- 
tive peculiar velocities inflow. Both redshifts and peculiar veloci- 
ties are in the CMB frame. Clusters are indicated by the CANs. 
Peculiar velocity errors are shown, but redshift errors (which arc 
small) are omitted for clarity. The dotted curves correspond to 
the typical itlcr peculiar velocity errors for clusters with peculiar 
velocities based on 3 galaxies. The unweighted mean peculiar ve- 
locity, and the number of sample clusters, are shown for the HCB 
and PPC regions separately and for the sample as a whole. The 
probability that the observed peculiar velocities are consistent 
with strictly zero motions is also given. 



liar velocity detection (J19, CAN=59) is omitted, this rises 
to 8%. If the peculiar velocity errors were under-estimated 
by 5% (10%), then the fit is consistent at the 6% (15%) 
level. If random thermal motions with an rms of 250kms~^ 
(500kms~^) are assumed, then the fit is consistent at the 
5% (30%) level. There is, therefore, no evidence in the EFAR 
sample for significant bulk motions in the HCB or PPC vol- 
umes. 

The components in Supergalactic coordinates of the 
mea n p eculiar velocity in redshift shells are shown in Fig- 
ure ha. There is no sign of any trend with redshift in the 
mean peculiar velocity, either for the whole sample or for 



the two regions separately. None of the components of the 
mean peculiar velocity are significant in any redshift bin 
apart from the 12000-13000 km s"^ bin in HCB, which is 
due to J19 (CAN=59)— cf. Figure |ll[ 

We can estimate the intrinsic dispersion of the pecu- 
liar velocity field using the maximum likelihood approach 
described in Paper VI (see Section 2.1 and Appendix A; 
cf. Watkins 1997). The upper panels of Figure [13 show the 
distributions of peculiar velocities, both radially and in Su- 
pergalactic coordinates, for the HCB and PPC regions sep- 
arately and for the whole sample. The peculiar velocities in 
all cases have means close to zero, and the question is how 
large an intrinsic dispersion is required, combined with the 
observational uncertainties, to reproduce the observed scat- 
ter in the peculiar velocities. The lower panels of Figure 
show the relative likelihood, Aln£ = ln£max — ln£, as a 
function of the assumed intrinsic dispersion. The most likely 
estimate of the three-dimensional velocity dispersion for the 
whole sample is about 600kms~^, but the Icr range is 0- 
1200 km s~^. The most likely dispersions for the HCB and 
PPC regions separately are about 300 km s~ ^ and 700 km s~ ^ 
respectively. Hence the intrinsic dispersion of the clusters' 
peculiar velocities is not well-determined by this data, due 
to the large uncertainties in the observed peculiar velocities. 



6.3 Comparisons with other results 

A comparison of the EFAR bulk motion to other measure- 
ments of bulk motions on various scales, and to theoretical 
predictions, is given in Figure The figure shows the re- 
ported bulk motions from a number of other observational 
studies as a function of the effective scale of the sample. 
Also shown is the theoretical prediction for the bulk motion 
measured with a top-hat window function of radius K (in 
h~^ Mpc) for a fairly 'standard' flat ACDM cosmology hav- 
ing a power spectrum with shape parameter r=0.25, nor- 
malisation (T8 = l-0 and Hubble constant h=0.7 (correspond- 
ing to ^0=0. 36 and f2A=0.64; see, e.g.. Coles & Lucchin 
1995, p.399). 

This comparison is limited by a number of factors: 
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Figure 13. Upper panels: The histograms of the pecuUar veloci- 
ties for the HCB and PPC regions and the whole sample, both ra- 
dially and projected in Supergalactic (X,Y,Z) coordinates. Lower 
panels: The relative likelihood, A In £ = In £max — In £, as a func- 
tion of the assumed intrinsic dispersion, both overall and in each 
Supergalactic coordinate. The solid curve is for the whole sam- 
ple; the dashed and long-dashed curves are for the HCB and PPC 
regions respectively. The upper dotted line is the \a confidence 
level for the whole sample, while the lower dotted line is the Icr 
confidence level for both individual regions. 

(i) The finite, sparse and non-uniform observed samples do 
not have top- hat window functions, and their effective scales 
R are not well-defined (compare this figure with the similar 
figure in Dekel (2000)); this uncertainty is ameliorated by 
the slow decrease in the expected bulk motion with scale. 

(ii) Only the amplitudes of the bulk motions are compared, 
and not the directions; however, the observed bulk motions 
that are significantly different from zero have a common di- 
rection to within about 30°, close to the direction of the 
CMB dipole. (iii) The uncertainties in the measured bulk 
motions are only crudely estimated in some studies, and ig- 
nore or under-estimate the systematic biases. Despite these 
limitations, the figure does show that, allowing for both ob- 
servational uncertainties and cosmic variance, the measured 
bulk motions are in most cases quite consistent with the 
theoretical predictions (which vary relatively little for any 
model that is consistent with the currently-accepted ranges 
of the cosmological parameters) . In this section and the next 
we determine the extent to which the EFAR results are con- 
sistent with the models and with the possibly-discrepant 
results of Lauer & Postman (1994; ACIF) and Hudson et al. 
(1999; SMAC). The bulk fiow obtained by Willick (1999; 
LPIOK) is similar to the SMAC result, and is not consid- 
ered explicitly. 

We can test whether the observed EFAR peculiar ve- 
locity field is consistent with the bulk motions claimed by 
other authors. The bulk motion of the Lauer & Postman 
(1994) cluster sample in the CMB frame, based on brightest 
cluster galaxy distances as re-analysed by Colless (1995), is 
764kms"^ in the direction (Z,&)=(341°,49°). This direction 
is only 39° from the median direction of the HCB clusters in 
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Figure 14. Bulk motion amplitude as a function of scale. The 
theoretical curve is the expectation for the bulk motion within a 
spherical volume of radius i? in a ACDM model (r=0.25, (T8=1.0, 
/i=0.7); the grey region shows the 90% range of cosmic scatter. 
The bulk motions determined in various studies are shown at the 
'effective scale' of each sample (which is generally only approx- 
imate). The bulk motions shown are for the Local Group w.r.t. 
the CMB (Kogut et al.l993), 7S (Lynden-Bell et al. 1988), ACIF 
(Lauer & Postman 1994; Colless 1995), SFI (Giovanelli et al. 
1998a), SCI (Giovanelli et al. 1998b), SCII (Dale et al. 1999a), 
Mklll (Dekel et al. 1999), SMAC (Hudson et al. 1999), LPIOK 
(Willick 1999), Shellflow (Courteau et al. 2000), SNe (Riess 2000), 
ENEAR (da Costa et al. 2000), SBF (Tonry et al. 2000), PT 
(Pierce & TuUy 2000), and EFAR (this work). Also shown are 
the predicted bulk motions derived from the PSCz redshift sur- 
vey (Saunders et al. 2000; Dekel 2000). 

the EFAR sample, and its antipole is just 15° from the me- 
dian direction of the PPC clusters. Consequently the EFAR 
sample is able to provide a strong test of the existence of 
the Lauer & Postman bulk motion. Figure [l^ shows the 
peculiar velocities of the EFAR sample as a function of the 
cosine of their angle with respect to the direction of the 
Lauer & Postman dipole. The best-fit bulk fiow in the Lauer 
& Postman direction has V=250±209 kms~^, and is consis- 
tent with zero at the 1.2(t level. A test finds that a pure 
Lauer & Postman bulk motion of 764 km s~^ in this direction 
is consistent with the data at only the 0.2% level. 

The bulk motion of the SMAC sample, for which pecu- 
liar velocities are derived from FP distances by Hudson et al. 
(1999), is 630 kms"^ in the direction (;,6)=(260°,-l°). The 
median direction of the HCB clusters is 57° from the an- 
tipole of this motion, and the median direction of the PPC 
clusters is 76° from the antipole. Hence the EFAR sample 
is less well-suited to testing for bulk motions in this direc- 
tion. Nonetheless, the formal rejection of the SMAC motion 
is even stronger than for the Lauer & Postman motion. Fig- 
ure shows the peculiar velocities of the EFAR sample 
as a function of the cosine of their angle with respect to 
the SMAC dipole. The best-fit bulk fiow along the SMAC 
direction has 1/=— 536±330 kms~^ (i.e. in the opposite di- 
rection), and is consistent with zero bulk motion at the 1.6(t 
level. A test finds that a pure SMAC bulk motion of 
630kms~^ in this direction is consistent with the data at 
only the 0.04% level. 

It is worth noting that an observed bulk flow amplitude 
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Figure 15. The peculiar velocities of the EFAR clusters versus 
the cosine of their angle with respect to the direction of (a) the 
Lauer & Postman dipole, {i,6)=(341° ,49°), and (b) the SMAC 
dipole, (i,fe)=(260°,-l°). Each cluster is indicated by its CAN. 
The solid line shows the claimed relation; the dotted line is the 
best fit to the EFAR data (see text for details). 



of zero would be consistent with the Lauer & Postman flow 
at less than the 0.2% level, but consistent with the SMAC 
flow at the 3.2% level — if the real bulk flow is small, there- 
fore, the apparently high significance of the rejection of the 
SMAC flow may be the result of the large uncertainty in the 
observed amplitude of the flow. 

These tests do not take into account the correlated 
errors in the peculiar velocity estimates. We therefore carry 
out Monte Carlo simulations of the EFAR dataset, includ- 
ing the effects of the correlated errors, in order to check 
the consistency of the observed peculiar velocities with the 
claimed bulk flows of Lauer & Postman (LP) and Hudson 
et al. (SMAC). Figure shows the distributions of the bulk 
flow amplitudes recovered from 500 simulations of the LP 
and SMAC bulk motions. The mean values of the recovered 
bulk flow amplitude (Vsim) are very close to the true values 
(Vlp or VsMAc), although in each case there is a small but 
statistically significant bias. However the value of the bulk 
flow amplitude derived from the actual EFAR dataset (Vobs) 
is in both cases far out on the wing of the distribution: only 
one of the 500 simulations of the Lauer & Postman flow, 
and none of the 500 simulations of the SMAC flow, yields a 
bulk flow amplitude less than the observed value. Hence the 
observations are consistent with a pure Lauer & Postman 
bulk flow only at the 0.2% level, and with a pure SMAC 
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Figure 16. Simulations of the recovery from the EFAR dataset 
of (a) the Lauer & Postman (1994; LP) bulk flow and (b) the 
Hudson et al. (1999; SMAC) bulk flow. The histograms are the 
distributions of the recovered bulk flow amplitude in the direc- 
tions of the LP and SMAC dipoles. The labelled arrows show the 
true amplitude (Vlp or Vsmac)i the mean of the recovered am- 
plitudes {Vsim) and its rms scatter, and the observed amplitude 
(Vobs) a'Ud its uncertainty. 



bulk flow at less than the 0.2% level. The correlated errors 
in the peculiar velocities do not significantly alter the results 



obtained from the x tests. 



6.4 Comparisons with theoretical models 

The above comparisons assume pure bulk flows and ignore 
the greater complexity of the real velocity field. We can 
make more realistic comparisons if we adopt a more detailed 
model for the velocity field. In principle this approach also 
allows us to use the observed peculiar velocities to discrim- 
inate between different cosmological models. The velocity 
field models are characterised by a mass power spectrum, 
which determines the velocity field on large scales where the 
dynamics are linear, and a small-scale rms 'thermal' motion, 
(T», which approximates the effects of non-linear dynamics 
on small scales. Given such a model, the method for com- 
puting the expected bulk flow in a particular sample, and 
for estimating the probability of an observed bulk flow, has 
been developed by Kaiser (1988) and Feldman & Watkins 
(1994, 1998). 

As shown by Feldman & Watkins (1994), the covari- 
ance matrix for the maximum likelihood estimator of the 
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bulk flow in a sample is given by the sum of a 'noise' term, 
which depends on the spatial distribution of the clusters, the 
errors in their peculiar velocities and the thermal rms mo- 
tions, and a 'velocity' term, which also depends on the power 
spectrum of the assumed cosmological model. We adopt a 
thermal rms motion of (t*=250 kms~^. Although this value 
is not well-determined it has little effect on the results (as 
we show below), since it enters in quadrature sum with the 
uncertainties on the cluster peculiar velocities, which are 
generally much larger (see Table ^). Our adopted cosmo- 
logical model has a CDM-like power spectrum with r=0.25 
and (78 = 1.0, consistent with the power spectrum measured 
from the APM galaxy survey (Baugh & Efstathiou 1993) 
and the PSCz redshift survey (Sutherland et al. 1999). This 
corresponds to the currently-favoured flat ACDM cosmology 
with Ho«70kms~^Mpc"\ r2o«0.35 and r2A~0.65. 

The survey's sensitivity to the power spectrum is deter- 
mined by its window function. Figure shows the window 
function for the EFAR sample along the Supergalactic X, Y 
and Z axes; the Y axis in particular shows the effect of corre- 
lated errors resulting from not having a full-sky sample. The 
model power spectrum is shown in Figure The product 
of the power spectrum and the window function, shown in 
Figure |l7p, gives the relative contributions of different scales 
to the covariance in the measured bulk velocity. The bulk 
velocity depends on a broad range of scales, with the largest 
contributions coming from scales of a few hundred h~^ Mpc. 

For the EFAR survey the 'noise' component of the co- 
variance matrix (in Supergalactic coordinates) is 



r [h"' Mpc] 
1000 100 



-fl01655 -f47914 -24001 
-f47914 +65373 -39617 
-24001 -39617 +87567 



while the 'velocity' component is 
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Figure 17. (a) The trace of the squared tensor window function 
for the EFAR sample along the Supergalactic X, Y and Z axes; 
(b) the power spectrum for a CDM-like model with r=0.25 and 
(78=1.0; and (c) the contributions of different scales to the covari- 
ance in the measured bulk velocity, given by the product of the 
power spectrum and the squared tensor window function. 



Thus the overall covariance matrix R is 



+138824 
+65125 
-24345 



+65125 
+88538 
-45701 



-24345 
-45701 
+108547 



.(15) 



It is immediately apparent that (for the model considered 
here) the covariance matrix is dominated by the 'noise' term. 

The maximum likelihood estimate, U , for the bulk flow 
of the sample clusters is given by 



(16) 



where Ui is the ith component of the bulk flow, Rlj is the 
'noise' covariance matrix, fnj is the jth component of the 
unit vector of the nth cluster, Vn and (t„ are the cluster's 
peculiar velocity and its uncertainty, and cr* is the assumed 
rms thermal motion of the model. For the EFAR sample, 
the maximum likelihood bulk flow vector in Supergalactic 
coordinates is (—24, —6, +717)kms~^, almost entirely in 
the SGZ axis. In Galactic coordinates this is 718kms~^ in 
the direction {l,b) = {45A° ,+5.9°). 

However this formal result is rather ill-determined, since 
it is far from the main axis of the EFAR sample (cf. Fig- 
ure ^). An indication of the uncertainty can be obtained 



by ignoring the cross-correlations in the covariance matrix 
and estimating the rms error as (Trace(i?'^))^/^=505 kms~^. 
In the context of the assumed cosmological model, the prob- 
ability of measuring a bulk flow vector U can be obtained 
by computing the statistic from the covariance matrix as 



(17) 



The probability (given the cosmological model and the prop- 
erties of the sample) of observing a bulk flow with a value of 
greater than this is given by the appropriate integral over 
the x^ distribution with 3 degrees of freedom (the 3 com- 
ponents of U). For the EFAR sample this procedure yields 
^■^=6.1 with 3 degrees of freedom, and hence the observed 
bulk flow is consistent with the model at the 11% confidence 
level. If the rms thermal motion a, is set to be zero rather 
than 250kms~^, the observations are still consistent with 
the model at the 9% confidence level. 

The expectation value for the bulk motion (given the 
cosmological model and the properties of the sample) can 
be obtained as 



V = 



((Ji(72(T3) 

(27r)»/2 



(18) 
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where ai, 02 and 173 are the lengths of the axes of the covari- 
ance eUipsoid obtained from the eigenvalues of the covari- 
ance matrix. The directions of these axes are given by the 
(orthogonal) eigenvectors of the covariance matrix. For the 
EFAR sample and our adopted cosmological model, these 
eigenvalues and eigenvectors (in Supergalactic coordinates) 
are 





= 454kms~^ 


ei 




= 309kms-i 
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0-3 


= 185 kms'^ 
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: (+0.7026, +0.5604, -0.4385) , 



(19) 



The corresponding directions in Galactic 
coordinates are ei = (172.6°,+30.6°), e2 = (82.0°,+1.2°), and 
e3=(350.0° ,+59.4°). We therefore find an expectation value 
for the amplitude of the bulk flow of 619 kms~^, so that the 
observed value is not much larger than that expected from 
our model, as the statistic indicates. It is worth noting 
that the expected bulk flow amplitude is strongly dominated 
by the 'noise' term in the covariance matrix. For our adopted 
cosmological model in the absence of noise, we would expect 
to measure a bulk flow amplitude from the EFAR sample of 
only 355 kms~^, whereas in the absence of any cosmological 
velocities, the noise in our measurement would still lead us 
to expect a bulk flow amplitude of 553kms~^. 

We obtain a smaller upper limit on the bulk motion 
we if consider only the component of the bulk flow along 
the minimum-variance axis of the covariance ellipsoid. Un- 
surprisingly, this axis, e3=(350.0° ,+59.4°), is just 20° away 
from the median axis of the 50 clusters in the peculiar ve- 
locity sample, ((/, 6)) = (7°,+42°). The expected bulk flow 
amplitude along this axis is 147kms-i (124 km s"^ from 
noise alone, 76kms~^ from model alone), while the max- 
imum likelihood estimate of the observed bulk motion is 
269kms~^. Since (73=185 kms~^, this gives x^=2.11 with 1 
degree of freedom, implying that the observed bulk motion 
in this direction is consistent with the model at the 15% 
confidence level. 

Thus there is no evidence that the bulk motion of the 
EFAR sample is inconsistent with a cosmological model hav- 
ing a CDM-like power spectrum with r=0.25 and (T8=1.0, 
consistent with the best current determinations. In fact, re- 
peating this analysis, we find that the observations are con- 
sistent with a wide range of cosmological models, including 
both standard CDM and open, low-density CDM models. 

We can also ask to what extent the EFAR sample is 
capable of testing whether the bulk motions measured by 
Lauer & Postman (1994), SMAC (Hudson et al. 1999) and 
LPIOK (Willick 1999) are consistent with the velocity field 
model. To do so we use the statistic computed according 
to equation inserting the EFAR covariance matrix for 
R and the observed Lauer & Postman, SMAC or LPIOK 
bulk motions for U . If the EFAR bulk motion had been 
found to be identical to the SMAC result, it would have been 
consistent with the velocity field model at the 25% level; if 
it had been found to be identical to the LPIOK result it 
would have been consistent with the model at the 9% level. 
However a bulk motion identical to the Lauer & Postman 
result would have been rejected at the 0.09% level. Hence, 
as expected, the directionality of the EFAR sample means 
that while it would have provided a strong indication of an 
inconsistency with the model if the Lauer & Postman result 




Figure 18. Contour plots of the bulk motion amplitude, in each 
direction on the sky, that would be rejected at the 1% level or 
better by the EFAR sample. The assumed power spectrum is 
CDM-like, with V = 0.25 and cr8=1.0, and the rms thermal mo- 
tions of the clusters is assumed to be tr, =250 km s~^. The EFAR 
clusters with positive (negative) peculiar velocities are indicated 
by asterisks (circles). Other symbols show the directions with re- 
spect to the CMB frame of the Local Group dipole (0), the Lauer 
& Postman (1994) dipole ((»), the SMAC (Hudson et al. 1999) 
dipole (®), and the LPIOK (Willick 1999) dipole (0). The con- 
tours run in steps of lOOkms"^ from 700kms~^ to 1500 km s~^, 
with the lowest contour being the thickest. 



had been recovered, recovery of the SMAC or LPIOK results 
would not have implied a problem with the model. 

We can generalise this analysis to illustrate how the di- 
rectionality of the EFAR sample affects the constraints it 
could place on observed bulk motions in different directions. 
Figure ^ shows, in each direction on the sky, the amplitude 
of the observed bulk motion that would be rejected as in- 
consistent with the velocity field model at the 1% confidence 
level using equation 

It is important to emphasise that although it would not 
have been surprising, under this model, to have recovered 
the SMAC motion from the EFAR sample, in fact the tests 
of the previous section indicated that the actual motions 
recovered from the EFAR sample are highly inconsistent 
with a pure SMAC bulk flow. As already noted, however, 
because those tests do not use a full velocity field model 
and do not account for the window function of the sample, 
they will tend to over-estimate the degree of inconsistency. 
The best test is a simultaneous consistency check between 
both datasets and the model (Watkins & Feldman 1995), 
determining the joint probability of deriving both the ob- 
served EFAR bulk motion from the EFAR sample and the 
observed SMAC motion from the SMAC sample under the 
assumptions of the velocity field model. This type of test has 
already been carried out for the SMAC sample with respect 
to various other samples by Hudson et al. (2000), who find 
consistency with all the other peculiar velocity surveys with 
the possible exception of Lauer & Postman, and a marginal 
confiict with a flat ACDM model similar to that used here. 
Once the SMAC peculiar velocities have been published, a 
similar test can be carried out to check the consistency of 
the EFAR and SMAC survey resuhs. 



7 CONCLUSIONS 

We have measured peculiar velocities for 84 clusters of galax- 
ies in two large, almost diametrically opposed, regions at 
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distances between 6000 and 15000 km s~^. These velocities 
are based on Fundamental Plane (FP) distance estimates for 
early-type galaxies in each cluster. We fit the FP to the best- 
studied 29 clusters using a maximum likelihood algorithm 
which takes account of both selection effects and measure- 
ment errors and yields FP parameters with smaller bias and 
variance than other fitting procedures. We obtain a best-fit 
FP with coefficients consistent with the best existing de- 
terminations. Apparent differences in the FPs obtained in 
previous studies can be reconciled by allowing for the biases 
imposed by the various fitting methods. We then fix the FP 
parameters at their best-fit values and derive distances for 
the whole cluster sample. The resulting peculiar velocities 
show no evidence for residual systematic errors, and, for the 
small numbers of clusters in common, are consistent with 
those measured by other authors. 

We have examined the bulk motion of the sample re- 
gions using the 50 clusters with the best-determined pe- 
culiar velocities. We find the bulk motions in both regions 
are small, and consistent with zero at about the 5% level. 
We use both direct comparison and the more sophis- 
ticated window function covariance analysis developed by 
Kaiser (1988) and Feldman & Watkins (1994, 1998) to com- 
pare our result with the predictions of standard cosmologi- 
cal models and the results of other studies. We find that the 
bulk motion of our sample is consistent (at about the 10% 
level) with the prediction of a ACDM model with parameters 
r=0.25, CT8=l-0 and h=0.1; indeed the motion is consistent 
with most cosmological models having parameters that are 
broadly consistent with the observed shape and normalisa- 
tion of the galaxy power spectrum. 

We examine whether our results can be reconciled with 
the large-amplitude bulk motions on similar scales found 
in some other studies. Our sample lies close to the direc- 
tion of the large-amplitude dipole motion claimed by Lauer 
& Postman (1994), so that we are able to make an effec- 
tive test of the bulk motion in this direction. We find that 
a pure Lauer & Postman bulk motion is inconsistent with 
our data at the 0.2% confidence level. This strong rejection 
of the Lauer & Postman result is supported by the win- 
dow function covariance analysis. We find an even stronger 
inconsistency between the EFAR peculiar velocities and the 
result of the SMAC survey (Hudson et al. 1999), with a pure 
SMAC bulk motion ruled out at the 0.04% confidence level. 
This is a surprisingly strong result, given that the main axis 
of the EFAR sample lies at a large angle to the direction 
of the SMAC dipole. It will be important to carry out a 
simultaneous consistency check of both datasets with a full 
velocity field model using the generalised covariance analysis 
described by Watkins & Feldman (1995) and Hudson et al. 
(2000). 

To summarise current observations of bulk motions on 
scales larger than 6000 kms~^ (i) The EFAR and SCII (Dale 
et al. 1999a) surveys find small bulk motions, close to the 
predictions of cosmological models that are constrained to 
be consistent with other large-scale structure observations. 

(ii) The SMAC survey (Hudson et al. 1999) finds a bulk 
motion with a much larger amplitude. However a full ac- 
counting for the uncertainties and window function of the 
survey shows that it is in fact only marginally inconsistent 
with the models (at about the 2a level; Hudson et al. 2000) . 

(iii) The LPIOK survey finds a bulk motion very similar to 



the SMAC dipole, but the smaller sample size means that 
the uncertainties are larger and consequently the result is 
not inconsistent, (iv) The Lauer & Postman (1994) result is 
inconsistent with such models at the 3-5% level (Feldman 
& Watkins 1994). However it is also inconsistent with the 
EFAR results (at the 0.2% confidence level) and with the 
other surveys combined (at the 0.6% level; Hudson et al. 
2000), and therefore should be treated with reserve. We con- 
clude that existing measurements of large-scale bulk motions 
provide no significant evidence against standard models for 
the formation of structure. 
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